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Abstract 

The Low Energy Telescopes on the Voyager spacecraft are used to 
measure the elemental composition (2 i Z i 26) and energy apectra (5 
to 15 MeV/nucleon) of solar energetic particles (SEPs) in seven large 
hare events. Four flare events are selected which have SEP abundance 
ratios approximately independent of energy/nucleon. The abundances 
for these events are compared from flare to fla^e and are compared to 
solar abundances from other sources: spectroscopy of the photosphere 
and corona, and solar wind measurements. 

The selected SEP composition results may be described by an 
average composition plus a systematic flare-to-flare deviation about 
the average. For each of the four events, the ratios of the SEP 
abundances to the four’-flare average SEP abundances are 

approximately monotonic functions of nuclear charge Z in the range 
6^ 2B. An exception to this Z-dependent trend occurs for He, whose 

abundance relative to Si is nearly the same in all four events. 

The four-flare average SEP composition is significantly difl'erent 
from the solar composition determined by photospheric 

spectroscopy: The elements C, N and 0 are depleted in SEPs by a factor 
of about five relative to the elements Na, Mg, Al, Si, Ca, Cr, Fe and Ni. 
For some elemental abundance ratios (eg. Mg/0), the difference 
between SEP and photospheric results is persistent from fiare to fiare 
and is. apparently not due to a systematic difference in SEP 
energy /nucleon spectra between the elements, nor to propagation 
effects which would result in a time— dependent abundance ratio in 
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individual flare events. 

The f cur-flare average SEP composition is in agreement with solar 
wind abundance results and with a number of recent coronal 
abundance measurements. The evidence for a common depletion of 
oxygen in SEPs, the corona and the solar wind relative to the 
photosph.'"s suggests that the SEPs originate in the corona and that 
both the SEPs and solar wind sample a coronal composition which is 
signiflcantly and persistently different from that of the photosphere. 
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Chapter 1 
Introduction 

The energetic particles which are often ejected from solar flares 
constitute a sample of solar material which may be analyzed to ]riel.d 
information on the elemental and isotopic composition of the solar 
atmosphere — information that impacts a wide range of 
astrophysical problems from the history of the solar system, to solar 
structure and dynamics, to nucleosynthesis. 

Following the pioneering spectroscopic study by H. N. Russell 
(1929), the composition of the sun has been the subject of extensive 
experimental and theoretical investigations (see. e.g.. Claas 1951, 
and Goldberg. Muller and Aller 1960), yet today our knowledge of solar 
composition remains inadequate. Currently, solar composition 
information is obtained using a number of techniques including 
spectroscopy of the photosphere and corona, and measurements of solar 
wind and of energetic particles from solar flares (see the review by 
Ross and Aller 1976. the update by Aller 1960 and references therein). 
None of the techniques are free of difllculty: Spectroscopic 
abundance determinations are (a) subject to inaccuracies in modeling 
the temperature and density structure of the solar atmosphere and 
to uncertainties in our knowledge of the spectral line formation 
mechanisms and atomic transition probabilities, and are (b) 
particularly diflicult to obtain for some elements — most notably 
helium, the second most abundant solar constituent. While the solar 
wind and energetic flare particles are samples of solar material whose 
composition is directly measured, there is the concern that these 


samplei maj be biased since in both cases the elemental corarpcsition 
varies and may be influenced by acceleration effects which are 
currently not understood. 


The potential of solar energetic particle (SEP) measurements as a 
source of solar composition information was first explored In the early 
1960's using rocket-borne nuclear emulsion experiments with energy 
thresholds near 40 MeV/nucleon. The early work (see the reviews 
by Biswas and Fichtel 1965; Bertsch, Fichtel and Reames 1972) suggested 
that, for nuclei with nearly the same nuclear charge to mass ratio 
(namely most of the abundant nuclei except protons), the SEPs might 
indeed represent an unbiased sample of solar material. In six different 
fiare events, the ratios of the numbers of Ke nuclei to medium group 
nuclei (M = C + N + 0), when counted in common energy/nucleon (i.e. 
particle speed) intervals, were found to be (a) approximately 
independent of the choice of the energy/nucleon interval and (b) nearly 
constant from flare to flare (Biswas, Fichtel and Guss 1962; Biswas et. 
al. 1963; Biswas, Fichtel and Guss 1966; Durgaprasad et. al. 
1968; Bertsch, Fichtel and Reames 1972). Combining the measured 
SEP He/M ratio (= 58± 5) with the spectroscopically determined M/H 
ratio, Bertsch, Fichtel and Reames (1972) obtained H/He = 16± 2, one 
of the few available estimates of the solar H/He abundance ratio. 
The early emulsion measurements achieved only poor statistical 
accuracy and charge resolution for the elements Ne, Mg, Si and Fe, but 
suggested that the SEP abundances of these elements relative to the M 
group nuclei were roughly constant from fl.are to flare and similar to 
spectroscopic abundance results (Biswas and Fichtel 1965; Bertsch, 


Fichiel and Reames 1972). In contrast to the apparent constancy 
of SEP composition for He and heavier nuclei, the SEP proton to He 
nuclei ratio was found to vary (by factors of 10 to 100) from flare to 
flare and as a function of time and energy/nucleon in single flare 
events (see, e.g.. Biswas and Fichtel 1865; Fichtel and McDonald 1967). 

The above results suggested that during acceleration and 
propagation to earth the SEPs were fully stripped of electrons and 
interacted mainly with large scale electromagnetic fields. In this case 
most of the abundant nuclei, except protons, would have nearly equal 
charge to mass ratios and would be affected almost identically by the 
electromagnetic fields, no matter how complex. However, the charge to 
mass ratio of protons would be a factor of two larger than that of the 
other abundant nuclei, and this difference might account for the 
observed variability of the SEP proton to He (or heavier) nuclei ratio. 
The early SEP results also suggested that processes which do not scale 
as the charge to mass ratio (such as thermal particle acceleration, 
or particle deceleration by Coulomb interactions) were not very 
important during SEP acceleration and propagation. 

However, following the early nuclear emulsion results, a large 
amount of research has indicated that, even for elemental species with 
the same nuclear charge to mass ratios, the SEP elemental composition 
does vary from flare to flare (see, e.g., Armstrong and Krimigis 1971; 
Armstrong et. al. 1972; Mogro-Campero and Simpson 1972b; Bertsch 
et. al. 1973; Teegarden, von Rosenvinge and McDonald 1973; Crawford 
et. al. 1975) and sometimes varies with time, space, and/or energy per 
nucleon in single flare events (see, e.g.. Bertsch, Biswas and Reames 
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1974: Van Allen. Venkeiarangan and Venkaieaan 1974; Crawford et. al. 
1975; Armstrong ei. al. 1976: O’Gallagher et. al. 1976; Scholer et. 
al. 1976). Some of the more extreme SEP elemental composition 
anomalies are associated with large ^He enhancements (Hovestadt et. 
al. 1975, Hurford et. el. 1975; Zwickl et. al. 1978). The few reported 
charge state measurements of SEPs are for energies below about 1 
MeV/nucleon and indicate that while C and 0 nuclei are almost fully 
stripped, iron nuclei retain more than half their electrons (Sciambi et. 
al. 1977; Gloeckler et. al. 1676). 

This apparent complexity, and the availabliity of high quality SEP 
composition measurements for only a small number of flare events, has 
made it diflicult to determine which features of SEP composition are the 
same from flare to flare and which are variable. One area of 
uncertainty is the energy dependence of SEP composition. Until 
recently, the most extensive work in this area was performed using 
plastic and glass track detectors and nuclear emulsions aboard rockets 
and Apollo spacecraft. Crawford et. al. (1975) summarized this work 
and concluded that, above an energy Eq (which ranged between about 5 
and 20 MeV/nucleon in five different flare events), the SEP composition 
was (a) approximately independent of energy/nucleon, (b) roughly 
constant from flare to flare and (c) the same as spectroscopic 
composition results within factors of two to three. However ' elow 
energy Eq, they found that the heavy elements were enha:.ce:.’ by an 
amount which varied from flare to flare, but always increased with 
decreasing energy/nucleon and increasing nuclear charge Z (see Figure 
1.1 which shows data from one of the five flares). In contrast, Mason, 


SEP composition measuremr:nts for the 25 January 1971 solar flare 
event, showing the enhancement of heavy elements, such as Fe, at 
energies below about 15 MeV/nucleon (taken from Crawford et. al. 
1975). The flux measurements of He nuclei by Lanzerotti, Maclennan 
and Graedel (1972) were performed with an instrument aboard the 
IMP 5 satellite. The other data, including those of Bertsch et. al. 
(1973), were obtained during a rocket flight. (The smooth curves 
drawn through the data are meant only to guide the eye.) 
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Hovestadt and Gloeckler (1979), in a satellite study of flare events in the 
1073 to 1977 period, found that the average SEP composition for the 
abundant elements from He through Fe was essentially the same at 1 
MeV/nucleon as that measured by others above 10 MeV/nucleon. 
Further, heavy element enhancements which were found only at low 
energies by Crawford et. al. (1975) have been reported above 25 
MeV/nucleon in one flare event (Bertsch and Reames 1977) and, in 
another event, were found to extend from 25 to above 100 MeV/nucleon, 
independent of energy/nucleon (Dietrich and Simpson 1976). 

Another area of uncertainty is the relationship of SEP 
composition to the solar composition values which have been determined 
by other means. Crawford ct, al. (1975) compared their best 
estimate of SEP composition above 15 MeV/nucleon with their best 
estimate of solar atmosphere composition based on spectroscopic data 
and concluded that the differences did not correlate with either first 
ionization potential or nuclear charge Z, and were small enough 
(factors of two to three) to result from errors in the spectroscopic 
abundances. On the other hand, Webber (1975) found that the 
differences between his best estimates of SEP and solar atmosphere 
composition were significant and were correlated with first ionization 
potential — the elements with high first ionization potentials were 
found to be depleted in SEPs (see Figure 1.2). Averaging over seven 
fiare events Mogro-Campero and Simpson (1972a,b) reported still 
different results — the SEP abundances of the heavy elements were 
enhanced relative to "solar system" composition (derived from a 
mixture of solar spectroscopic data, meteoritic data, earlier SEP 


F'^wre 1J3 


A comparison of solar energetic particle (SEP) abundances to solar 
atmosphere (SA) abundances determined by spectroscopy, showing a 
relative depletion in SEPs of elements with high first ionization 
potential (taken from Webber 1975). 
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results, and nucleosynthesis theory: Cameron 1968) by an amount 
which increased as an approximately monotonic function of nuclear 
charge Z, for the abundant nuclei from carbon through iron. 
Similar heavy element enhancements were reported by Bertsch and 
Reames (1677) and Dietrich and Simpson (197B); see Figure 1.3. 
Furthermore, Dietrich and Simpson (1978) claimed an additional 
enhancement of the rare odd-Z nuclei, such as B, F, Na. and Al, by 
amounts consistent with the production of these nuclei in spallation 
reactions of the other even-Z nuclei during passage through an 
estimated 0.6 gm/cro® of solar atmospheric material — a claim which 
has been disputed by Cook et. al. (1976) and McGuire, von Rosenvinge 
and McDonald (1679) (also see Section 4.4.2 of this thesis). 

The above comments indicate the need to study SEP composition 
in a large number of flare events with a single, high-quality 
experiment and common analysis criteria, as is done in this thesis. 
Here, we present SEP composition measurements for seven large flare 
events which occurred in the September 1977 to May 1970 period. The 
measurements were performed in interplanetary space with advanced 
instrumentation on the Voyager 1 and 2 spacecraft and cover the 
nuclear charge range from He (Z=2) through Ni (Z=28) with the best 
combination of charge resolution, background rejection and statistical 
accuracy achieved to date. Thus, for the first time, we are able to study 
statistically significant measurements of 15 different elemental species 
(He, C. N, 0, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Fe, and Ni) in seven individual 
flare events. 


With the goal of understanding the relationship between SEP 


Figure 


The ratio (Enhancement Factor, Q) of solar energetic particle (SEP) 
abundance to 'holar system” abundance plotted versus nuclear charge Z, 
for the 24 September 1977 flare event (taken from Dietrich and Simpson 
1978). The use of "solar system” abundances which were a mixture of 
results from difl^erent sources (including earlier SEP data for the 
elements Ne, Mg, and Si) led Dietrich and Simpson to claim that (a) the 
enhancement of even-Z nuclei in SEPs increased as an 
approximately monotonic function of Z and (b) the rarer odd-Z nuclei 
Na and A1 were additionally enhanced in SEPs as the result of 
spallation of even-Z nuclei in the solar atmosphere. Cook et. al. (1979) 
and McGuire, von Rosenvinge and McDonald (1979) obtained difl'erent 
conclusions (see section 4.4.2 of this thesis). 
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elemental composition and the composition of the sun we proceed in 
this thesis according to the following scheme; 

(1) We select only large flare events. In addition to affording 
the necessary statistical accuracy, large flare events show less SEP 
compositional variability than smaller events (Z^rickl et. al. 1976; Mason, 
Hovestadt and Gloeckler 1979). Furthermore, extreme SEP 
compositional anomalies associated with ^He enhancements appear to 
occur only in relatively email events (Hurford et. al. 1975. Zwickl et. 
al. 1976). Thus, while the study of smaller events should ultimately 
provide insight into details of flare acceleration mechanisms and 
small scale solar atmospheric inhomogeneities, the larger events 
(where energetic particle acceleration may occur over large portions of 
the solar atmosphere) are more suitable for the study of global 
composition. 

(2) We attempt to minimize the effects of possible SEP 
"acceleration/propagation bias" by (a) the selection of averaging time 
periods for each flare which exclude times when the composition is 
likely to be affected by propagation and (b) the rejection of those flare 
events in which the measured SEP composition is dependent on 
energy/nucleon. Here, the term "acceleration/propagation bias" refers 
to any effect which causes the SEP elemental composition, in a fixed 
energy /nucleon interval, to differ from the composition of the 
pre-accelerated plasma at the SEP fiare acceleration site. By 
examining the measured SEP composition as a function of time and 
energy/nucleon in each fiare event we seek to both gauge the extent of 
possible acceleration/ propagation biases and minimize their effect. 
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(3) We compare SEP composition measurements among the 
different flare events. 

(4) We compare our SEP composition measurements to solar 
abundance measurements from other sources — (a) spectroscopic 
results from the photosphere, (b) spectroscopic results for the corona, 
and (c) solar wind measurements. 

The main result is the discovery of large persistent 
differences between SEP elemental composition and the results of 
photospheric spectroscopy (which are currently considered thf? most 
reliable of the different types of solar abundance measurements). The 
differences are apparently not due to any acceleration/propagation bias 
which we could detect in our energy range, but are roughly ordered 
by first ionization potential — an atomic parameter which plays a 
major role in the physics of the photosphere. These results suggest that 
(a) the composition of the SEP acceleration site is persistently 
different from that measured for the photosphere and/or (b) there is a 
persistent SEP acceleration bias operating at energies well outside our 
range of observation. The agreement we find among our SEP elemental 
composition results, solar wind data and a number of recent coronal 
abundance measurements favors possibility (a) and suggests 
(assuming the correctness of the photospheric composition results) 
that the SEPs originate in the corona, sampling a coronal composition 
which is significantly different from that of the photosphere. 

Brief accounts of this work have been published in the 
proceedings of the International Cosmic Ray Conference, Kyoto, Japan 
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(Cook et. al. 1979) and in the Astrophysical Journal (Letters) (Cook. 
Stone and Vogt, 1960). In addition, Mewaldt (1980) reviewed recent 
progress in SEP elemental and isotopic composition 

measurements, including preliminary results of this work. 
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Chapter 8 
The Experiment 


2.1 The Spacecraft 

The observations reported in this thesis were performed with the 
Low Energy Telescopes (LETs) of the cosmic ray detector systems 
aboard the Voyager 1 and 2 spacecraft. The Voyagers were launched 
toward the outer heliosphere in the fall of 1977 and followed 
trajectories nearly in the ecliptic plane (see Figure 2.1). The 
spacecraft are three-axis stabilized. All of the observations 
reported here were obtained during the September 1977 to May 1978 
perioa when both spa<;e'^raft were in interplanetary space, well beyond 
the influence of the earth's magnetic field. 

2.2 The Voyager Cosmic Ray Subsysteme 

The essentially identical Cosmic Ray Subsystems (CRS, see Stone 
et. al. 1977) on Voyagers 1 and 2 each consist of four Low Energy 
Telescopes (LETs), two High Energy Telescopes (HETs), the Electron 
Telescope (TET), and their associated electronic data systems. Together 
these telescopes measure the energy spectrum of electrons from 3-110 
MeV and the energy spectra and elemental composition of nuclei from 
hydrogen through nickel over an energy range from 3-500 
MeV/nucleon. The telescopes, shown in schematic cross section in 
Figure 2.2, employ silicon solid-state detectors exciusively and are 
designed to achieve (1) the reliability required for the anticipated 20 
year mission life, (2) exceUent charge and energy resolution, and (3) 


Pifure 2.1 


The trajectories of the Voyager 1 and 2 spacecraft from launch to 
Jupiter encounter. The data for this work were obtained in the 
September 1977 to May 1078 time period. 
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Figure 2J3 

Schematic cross sections of the particle telescopes of the Voyager 
Cosmic Ray Subsystem. 


r 
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the high count rate capability necessary during large solar flares 
and passage through the Jovian and Saturnian magnetospheres. 

The Voyager CRS is the result of a collaboration of scientists and 
engineers at the California Institute of Technology, the Goddard 
Space Flight Center, the University of Arizona, and the University of New 
Hampshire. 

2.3 The Low Energy Telescope System 

The LET system on each Voyager spacecraft incorporates four 

nominally identical charged particle telescopes (A, B, C and D) which 

ei F 

use the -E technique (described in the next section) to measure the 

kinetic energy and the nuclear charge Z of individual Incident nuclei in 
the range 1^2^26. The kinetic energy range of response varies 
from about 3-0 MeV/nucleon for protons and helium nuclei to about 
5-30 MeV/nucieon for iron nuclei. The four telescopes, in addition to 
having a relatively large combined geometry factor of 1.7 
cm^ steradian, are oriented at different viewing angles to provide 
three dimensional information on energetic particle streaming 
patterns. 

Each LET (see Figure 2.2) contains four totally depleted silicon 
surface barrier detectors, labeled L1-L4. LI and L2 are nominally 
identical 35 fjsn. thick detectors of the 'keyhole" design, where the active 
area is precisely defined by the location of vapor-deposited aluminum 
and gold contacts of about two centimeters diameter. L3 and L4 are 450 
/im thick detectors with active areas about 2.4 centimeters in 
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diameter. A thin (3//m) aluminum foil covers the top of each LET, 
protecting the detectors from sunlight and providing thermal control. 

The detectors L1-L3 are connected through charge sensitive 
preamplifiers and shaping amplifiers to linear 4096 channel pulse height 
analyzers (see Figure 2.3). Threshold circuits connected to the L1-L4 
amplifier outputs provide digital signals to coincidence circuitry which 
determines when pulse height analysis occurs. Pulse height analyzed 
(PHA) events are automatically sorted into two groups, Z< 3 (protons 
and alphas) and 3 (lithium through nickel), by the discrimination of 
an appropriate linear combination of the LI, L2 and L3 analog 
signals. The PHA events are temporarily stored in separate Z< 3 and 
Z^ 3 buffers which are sequentially polled for readout into the Voyager 
telemetry stream. Thus, during times of high counting rate (e.g. major 
solar flare events) when the PHA event readout rate is limited by 
telemetry, the occurrence of the relatively rare ZS 3 PHA events is 
enhanced in a predictable way. A rate accumulator system monitors 
single detector and various coincidence counting rates, including 
those needed to normalize the PHA event sample to obtain absolute 
flux measurements. (Details on the CRS electronic data system are 
given in Stilwell et. al. 1979) 

Throughout the observation periods reported here the basic 
requirement for pulse height analysis (which can be controHed by 
ground command) was a coincidence of discriminator signals from 
detectors LI and L2. L4 was at all times in anticoincidence so that 
normally only particles entering LI and stopping in either L2 or L3 
were pulse he.ght analyzed. 
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Figure 2.3 

Schematic diagram of one HET/LET electronic system. An identical 
system serves the other HET and two LETs. Note that the HET and LET 
systems share post amplifiers and pulse height analyzers. 
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The nominal characteristics of detectors L1-L4, and their 
corresponding discriminator thresholds and PHA channel widths are 
listed in Table 2.1. (The actual values of most of these parameters 
were measured as discussed in the calibration section of this chapter.) 

2.4 The -4~ ~E Technique 
dx 

The charge and energy of individual incident nuclei are measured 
ciE 

using the technique as illustrated in Figure 2.4. An incident 

nucleus (nuclear charge Z, mass M, and velocity v) penetrates a thin 
front detector of thickness L, and then deposits most of its initial 
energy and stops in a second detector. The (non -relativistic) 
ionization energy loss AE measured in the thin detector is roughly 
proportional to LZ^/v^ while the total energy E=Mv^/2 is the sum of AE 
and the residual energy E' measured in the second detector. Therefore, 
the product AExE is roughly proportional to LZ^M and is strongly 
dependent on Z. The response of a typical LET is illustrated in Figure 
2.5, a plot of the energy deposited in LI (AE) versus the sum of the 
energies deposited in L2 and L3 (E') for a raw sample of Z^ 3 events from 
LET C of Voyager 1. The "tracks" of the relatively abundant elements C, 
N, 0, Ne. Mg, Si, S and Fe are apparent, as are the less populated tracks 
of Na, Al, Ar and Ca. The finite width of the tracks is mainly due to the 
variation of particle incidence angles and detector thickness 
non-uniformities, which both contribute to a variation of the 
pathlength, L. 


For nuclei which penetrate both detectors LI and L2, and then stop 
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^-E TECHNIQUE 



Approximate Calculation: 

AE oc LzVv^ 1 ^ AE-E a LZ^M 

AE + E' = E “ J 

In practice, the charge, Z, is calculated by solving: 


L = R(AE + E'’,Z,M)-R(E'Z,M) 
with M = f(Z) a 2Z 


Figure 2.4 


2B 


Fi^re 2.5 

dE 

The — -E response of a t 3 rpicai LET. The energy deposited in detector LI 

(AE) is plotted versus the sum of the energies deposited in detectors L2 
and L3 (E'). for a sample of 3 PHA events from LET C of Voyager 1. To 
prevent plot saturation, only every tenth event was plotted for elements 
oxygen and below. 



500 
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in L3, two independent measurements are obtained. This redundant 
information improves the charge resolution and significantly reduces 
background for the rarer elements (as will be seen in Section 3.5). 

2.5 Calibrations 

2.5.1 Detector Thicknesses and Areas 

The thickness and area of each LET detector was measured as 
described in detail by Gehrels (1980). The measurements are briefly 
discussed below, 

The thickness (« 35 /xm) of each of the Ll and L2 detectors was 
determined using laboratory measurements of the energy deposited by 
penetrating 0.70 MeV alpha particles from a Po source, together 
with the alpha particle range-energy relation of Vidor (1975). For 
each detector, a sequence of aluminum masks was used to obtain 
separate exposures of six concentric, but non-overlapping, annular 
regions, which together covered the entire detector. The mean of the 
alpha particle energy loss distribution from each annular region was 
used to determine the average thickness of the annulus, while the 
spread of the distribution was used to estimate the rms of the thickness 
variations of the annular region. These data were then used to 
compute the mean pathlength in the detector and the rms variation 
in the pathlength expected for particles of an isotropic flux 
penetrating both Ll and L2. The mean pathlength measurements 
were reproducible with a .'sigma of about 0.2 percent. The main 
systematic error is due to uncertainty in the alpha particle 
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range-energy relation. For example, use of the range-energy 
relation obtained by scaling the proton range-energy relation of Janni 
(1966) gives 3 percent larger mean pathlengths than those listed in 
Gehrels (1980), while the alpha range-energy relation of Ziegler (1977) 
gives 6 percent larger results. 

The thickness (« 450 tim) of each of the L3 and L4 detectors was 
measured relative to a standard gauge-block using capacitive probes 
and a precision micrometer. For each detector the thickness was 
measured at the center and at four positions around the edge, and the 
mean pathiength and rms variation of the pathlength for an isotropic 
flux of particles penetrating both Ll and L2 was calculated by 
assuming that the detector surfaces were, to first approximation, 
spherical sections. The thickness measurements were reproducible 
with about a 2 /^m sigma. 

A travelling microscope was used to measure the dimensions of the 
A1 and Au contacts which define the active areas of the Ll and L2 
detectors. In addition, the relative areas of all Ll and L2 detectors were 
obtained more precisely by placing each detector a standard 
distance from an ^*^^Am source and measuring the counting rate of 5.5 
MeV alpha particles. These data and the L1-L2 separation distances 
measured during the LET assembly process were used to calculate 
the LET geometry factors listed in Gehrels (1980). 

2.5.2 Detector Dead Layers 


For silicon surface barrier detectors the thickness of the dead 
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layers is very nearly equal to the thickness of the front (Au) and back 
(Al) electrodes, These thicknesses were specified for the LET detectors 
by their manufacturer, ORTEC, and were all within the range 40. 0± 0.6 
/xg/cm®. Since this is only about 0,5 percent of the total thickness of 
the LET 35 ^J.m detectors, the effect of the dead layers was not 
explicitly included in the particle energy measurements discussed 
later. 

2.5.3 LET Window Thicknesses 

The thickness of the thin Al window at the top of each LET was not 
measured in the laboratory. However, in determining particle 

incidence energies, a correction was made for energy loss in a window of 
thickness 3 the thickness specified by the manufacturer. The 

energy correction was largest (about 8 percent) for nickel nuclei near 
their two detector threshold of 5 MeV/nucleon, but deci ;ased rapidly 
with increasing energy and was only 3.5 percent for nickel nuclei at 
their three detector threshold of 8.7 MeV/nucleon, 

2.5.4 Electronic Energy Calibration 

The LET preamp-postamp-pulse height analyzer chains were 
calibrated in the laboratory as described in Povlis (1980). A 
precision pulser was used to inject a charge pulse into a preamp input 
and the corresponding pulse height was read out. For each LET 
detector i, the input charge amplitude Qj corresponding to a channel 
threshold P was obtained for about 20 selected channels at two 
temperatures; T = 0 °C and 20 °C For each detector the charge 
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Qi(P ,T) was linear in P to within 2 percent of full scale and Q|(P ,0) was 1 
to 2 percent larger than Qi(P,20). An estimate Ecay of the energy 
which, when deposited in LET detector i by an incident nucleus, would 
give a pulse height, P, was obtained using; 

Ecai.i(P.T)*E,Q,(P.T) (2.1) 

where E, is ratio of the ionization energy deposited to charge output 
for silicon detectors. The main uncertainty in this calibration was 
in the absolute amplitude of the calibrating charge pulses^; the 
uncertainty in their relative amplitudes was negligible due to the 
linearity and zero offset of the precision pulser. Also negligible, in 
the LETs’ energy range, is the error which results from the 
approximation that the charge output of silicon detectors is 
proportional to the energy deposited . Thus, the relation between 
energy deposited and pulse height was taken to be: 


Ei(P,T) = FrEc,u(P.T) 

wh'Te tl 0 Fj were determined using oxygv. ?HA events as described 
in the next chapter. 


1. The uncertainty in the absolute amplitude of the calibration 
charge pulses was due to the uncertainty in the capacitance of the 
test capacitors (a different one for each preamp) used to couple the 
precision pulser to the preamp inputs 

2. The pulse height response characteristics for heavy ions in silicon 
surface barrier detectors hes been studied (Kaufman el. al. 1974, and 
references therein) by comparing the true energy E|_ of incident heavy 
ions of carbon through uranium to the alpha particle energy Ea yielding 
the seme pulse height. The energy defect, Ej-Ea, was essentially zero for 
carbon and oxygen ions but increased with nuclear charge and was 
about 2 MeV for nickel ions near 1 MeV/nucleon. The size of the defect 
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was roughly independent of energy, suggesting that at the LETs* 
threshold for nickel (about 5 WcV/nucleon) the percentage defect 
would be only about 1 percent. 
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Chapter 3 
Data Analysis 


3.1 Orerriew 

The various steps in the data analysis are shown schematically in 
Figure 3.1 and include (1) the completion of the energy versus pulse 
height calibrations using oxygen PHA events acquired in-flight, (2) the 
use of in-flight neon, magnesium, silicon and iron PHA events to obtain 
a charge calibration of the experiment, (3) the selection of PHA events 

dE 

by 'two parameter" (-^ E) analysis for flux and abundance 

dx 

measurements of the elements C, 0, Ne, Mg, Si, and Fe over the full 

LET energy range, and (4) the selection of PHA events by 'three 

H P H F* 

parameter" (-^ dlc~^^ analysis for very low background flux and 

abundance measurements of helium through nickel. 

3.2 In-flight Energy Calibration 

As discussed earlier, the relation between energy, Ej and pulse 
height, P. for each LET detector, i. was taken to be: 


E,(P.T) =F,E„,j(P.T). (3.1) 

The functions Ecai.i were measured (as discussed in Section 2.5.4) at 
about twenty selected channel thresholds and two different 
temperatures during the laboratory electronic calibration and, 
throughout this work, were extended to intermediate valu»^.s of P and T 
by linear interpolation. The Fj were obtained by least squares fits 
which optimized the agreement between the locations of the in-flight 
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oxygen "tracks" (e.g.. see Figure 2.5) and those calculated using 
equation (3.1) together with the oxygen range-energy relation of Vidor 
(1975) and the LI and L2 detector mean pathlength measurements 
discussed earlier. The method is described below, using LET C of Voyager 
1 as an example. 

Figure 3.2 illustrates the iterative technique used to obtain the 
location of the oxygen track on a plot of L2 versus L3 pulse height. In 
the first pass (Figure 3.2a), candidate oxygen PHA events were selected 
using curves which loosely bracketed the track. The selected events 
viere used to define a preliminary oxygen track location in tabular 
form (Figure 3.2b) by calculating the mean pulse heights (Pl2 . Pl3)j and 
the standard deviation of the L2 pulss heights Sj, for each group j, of 
Nj (*« 20) PHA events with adjacent L3 pulse heights. The Sj were 
smoothed by averaging each consecutive set of ten values along the 
track, and the smoothed values were used to define the final selection 
curves at about ± 2.5 sigma from the tabulated track center. The 
events selected using these curves are shown in Figure 3.2c and were 
used to calculate the final track shown in Figure 3. 2d. The Ll versus L3 
oxygen track was obtained in a similar way. 

The first step of fitting the L2 versus L3 oxygen track was to 
convert the mean pulse heights and sigmas to energies using the 
laboratory electronics calibration, i.e. equation (2.1): 

(Pl2 • Pl 3 . S)j -♦ (Ecal.LS - Ecal,L3 • ^ )\ 

where Cj is the statistical uncertainty in Eeai,L 2 ,j- Then, for oxygen 
nuclei, the relationship between the average energy Et^, deposited in L2 
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Figure 3.2 

The selection of PHA events from flight data to define the location of the 
oxygen track (see text). 
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and the average energy Ei^s. deposited in L3 was calculated using; 

Ei2 — R”*[R(El 3) Tug] 

» f(^L3) 

where R is the oxygen range-enc»gy relation of Vidor (1975), and T^g is 
the L2 detector mean pathlength (see Section 2.5.1). Finally, F^g and Fi;j 
could be found by minimizing: 

X.«.(Fu.Fu) • 

J 

Similarly, F^j and a second determination of Fu could be obtained from 
the LI versus L3 oxygen track. 

For the Voyager 2 LETs, the Fj determined for two different flare 
periods (September 19 through 27 of 1977 and February 13 through 20 
of 1978) were found to differ by less than 1 percent. However, the 
values of Fls from the Ll versus L3 tracks were systematically 
larger then those obtained from the L2 versus L3 tracks by about 2 
percent. The finally adopted values of Fj (listed in Table 3.1) were 
obtained by minimizing: 

x“(Ft,, , F.^ . Fu) = [xf,(Fu . F^a) x!,(Ful ■ Fas)] 

period 1 

+ [xf3(F"Ll ' + Xa^(F’L2 • ^'u)] 

periods 

The good agreement obtained between the oxygen tracks and the fits 
calculated using the adopted Fj is shown in Figure 3.3, for LET C of 
Voyager ’ . 
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TABLE 3.1 

LET Gain Correction Factors 


^L1 ^12 ^L3 • 

Voyager 1 


LET A 

0.9719 

0.9864 

0.9278 

LET B 

0.9893 

0.9921 

0.8930 

LET C 

1.0320 

1.0164 

0.9035 

LET 0 

1.0010 

0.9733 

0.9695 


Voyager 2 


LET A 

0.9421 

0.9740 

0.9303 

LET B 

0.9836 

0.9561 

0.9007 

LET C 

0.9519 

1.0124 

0.9340 

LET D 

1.0049 

1.0406 

0.9377 
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Fifure 3.3 

Comparison of the oxygen track obtained from flight data (DATA) and 
the calculated oxygen track (FIT) for LET C of Voyager 1 (sec text). 
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Since the energj calibration of the L3 detectors was performed 
without making use of the L3 detector thicknesses (which were 
accurately measured in the laboratory; see Section 2.5.1) it was 
possible to obtain an independent check on the possible systematic 
errors accumulated during the entire energy calibration procedure. This 
was done as illustrated in Figure 3.4 which shows a histogram of the 
energy deposited by oxygen nuclei in detector L3 of Voyager 1 LET C. 
The high energy cutoff results from the penetration of oxygen nuclei into 
the anticoincidence detector L4. For each L3 detector, the predicted 
cutoff energy Eg (calculated using the oxygen range-energy relation of 
Vidor 1975 together with the L3 detector mean pathlength 
measurements) was compared to the observed cutoff energy energy Eg 
(the energy at which the response has roughly halved). The observed 
cutoff energies, Eg , were systematically about 3 percent lower than 
the predicted cutoff energies, , indicating that absolute systematic 
errors on the order of 3 percent were accumulated in the calibration 
procedure. The effect of such possible errors on the observations is 
discussed in Chapter 4. 

The rrns uncertainties due the finite number of PHA events used to 
define the oxygen tracks are only about 0.1 percent for the F^j and F^j 
and 0.2 percent for the 

3.3 The Charge Calibration 

For each PHA event, the energies .E^j, Eiz and E^ (if L3 was 
triggered) were obtained from the corresponding LI, L2 and L3 detector 
pulse heights using the energy calibrations established in the 
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Figure 3.4 

A histogram of the energy deposited by oxygen nuclei in detector L3 of 
Voyager 1 LET C. The histogram cuts off at high energies due to the 
penetration of oxygen nuclei through L3 into anti-coincidence detector 
L4. Ec and Eg are the observed and predicted cutoff energies 
respectively (see text). 
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previous section. For every event, a charge measurement Z1 was 
calculated by numerically solving; 

TlI * . Zl , M) “ R(Ejjsi+El 3 . Zl , M) 

where is the mean pathlength for the appropriate LI detector, 
R(E,Z,M) is the range in silicon of a nucleus with energy E, nuclear 
charge Zl, and mass M. For the purpose of solving the above equation, 
the mass, M, was approximated by the following continuous function of 
Zl: 


r 2(Z1) for Zl g 20.0 

M s \ 40.0 + 4.772(Zl-20.0) for 20.0 g Zl < 21.0 
( 2.132(Z1) for Zl ^ 21.0 

(This approximation is accurate for most of the abundant isotopes from 
carbon through iron — the small errors which are incurred for some of 
the rarer isotopes produce predictable shifts in the calculated charge 
measurement which are negligible for this study.) 

Similarly, for three detector events (where the L3 detector was 
triggered) a second charge measurement Z2 was calculated by 
solving: 


Ti^ = R(El2+El3 . Z2 . M) - R(El3 . Z2 , M). 

The generalized range-energy function, R, was taken to be of the 
semi-empirical form used by Heckman et. al. (1960) to summarize 
their measurements of the range of heavy ions in nuclear emulsion; 


48 


R(1!.M.Z). (i)Rp{|) . 

The first term in the expression is the particle range as scaled 
from the proton range-energy relation, Rp. of Janni ( 1966). The second 
term corrects for charge-pickup which becomes significant in the 
LETs energy range for high Z nuclei such as iron. The function C, 
appropriate for silicon detectors, was obtained from the Voyager 2 
flight data as follows; 

(1) The 'Ll versus L2+’L3" and 'L2 versus L3" energy loss tracks 
of neon, magnesium, silicon, and iron were derived from the Voyager 2 
flight data using an iterative technique like that described earlier for the 
oxygen tracks. 

(2) Guided by the results of Heckman et. al. (1960) the function 
C was parameterized as: 

C(x)=A,A2[exp(-.i-^)-l.o] 

+ A, [a 4 - Aglnfl.O + I . 

(3) The parameters Ai were adjusted using a computer program to 
get the best agreement between the Ne, Mg, Si and Fe tracks obtained 
from the flight data and the corresponding track locations calculated 
using the above generalized range-energy relation and the Ll and L2 
detector mean pathlengthc. The resulting fits were not perfect and 
it was necessary to weight the Fe data more heavily to obtain an 
adequate fit for Fe at the expense of Ne, Mg and Si. The adopted 
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function C is plotted in Figure 3.5. 

The results of the charge calculations are illustrated in Figure 
3.6, a scatter plot of the charge measurement Z1 versus total energy, 
Elj+El 2 +El 3 . for a typical LET. The calculated charge is adequately 
independent of energy, except near threshold. A similar result was seen 
for the other LETs and on plots of Z2 versus energy. (The unusual 
clump of events seen in Figure 3.6 at the energy threshold near Zl=2B is 
due to a background effect which will be discussed later.) 

3.4 Two Parameter ( 4 ^ -E) Analysis 

dx 

The selection of PHA events for use in the 'two parameter" 
abundance and flux measurements for nuclei of C, 0, Ne, Mg, Si and Fe 
(presented in Section 4.4) was made as follows: 

(1) The events were sorted according to their Zl charge 
measurement using charge boundaries (listed in Table 3.2) chosen to 
lie at about ± 2.5 sigma from the mean value of Zl for each 
element, thus excluding a negligible fraction of "good" events. 

(2) For each element, the events were then sorted by their total 
energy measurement, Eli+El 2 +El 3 , into bins corresponding to incident 
kinetic energy/nucleon intervals (listed in Table 4.4) that were chosen 
to be within the energy response range of all the LETs used. The total 
energy (Eli+Ei^+Eu) interval corresponding to a given incident 
kinetic energy/nucleon interval was calculated using the mean nuclear 
mass for each element (averaged over the isotope fractions given by 
Cameron 1973) and a small correction for ionization energy loss in the 


50 


Figure 3.5 

The adopted range— energy extension function, C(x). (The values of the 
parameters A1 through A6 which define this function are 3.2209, 2.0040, 
0.92383, 1.1242, 0.28267, 0.98002, respectively — see text.) 
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Figure 3.6 

A scatter plot of the charge measurement Zl versus energy (Eh+Ei^j+Els) 
for PHA events obtained with LET B of Voyager 2 during the period 
September 1977 to May 1978. (To prevent plot saturation, only every 
tenth event is plotted for Zl < 9. PHA events with < 2 MeV or Zl < 4 
are primarily due to background effects that will be discussed in 
section 3 , 5.2 and are not included in this plot.) 
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TABLE 3.2 

Charge Boundaries for 
Two Parameter (Zl) Analysis 


z 

Element 

Lower 

Boundary 

Uoper 

Boundary 

6 

C 

5.78 

6.27 

8 

0 

7.67 

8.36 

10 

Ne 

9.63 

10.42 

12 

Mg 

11.49 

12.50 

14 

Si 

13.35 

14.53 

26 

Fe 

(see Appendix- B) 

- 
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thin A1 entrance windows. The term '’energy/nucleon", as used here 
and throughout this work, actually refers to energy per proton mass 
unit. Since, for nuclei of H through Ni the difference between the 
number of nucleons and the mass in proton mass units is less than 1 
percent, this error in terminology in negligible. 

The two parameter analysis of iron nuclei required special 
attention (described in Appendix B) due to the background mechanism 
which caused the unusual clump of events seen in Figure 3.6 near 
Zl=28. 

, dF HF 

3.5 Three Parameter (■§ E) Analysis for Z§ 3 Nuclei 

dx dx 

The selection of PHA events for use in the low background 
abundance measurements for nuclei of lithium through nickel 
(presented in Section 4.3) was made as follows: 

(1) For each PHA event the charge measurements Zl and Z2 were 
checked for consistency as described in Section 3.5.1 below. 

(2) PHA events with consistent charge measurements were sorted 
into element bins using the average of Zl and Z2, and the charge 
boundaries discussed in Section 3.5.1 below (and listed in Table 3.3). 

(3) PHA events in each of the element bins from carbon through 
nickel were then counted if the total energy measurements, 
Eli + El 2 +El 3 , were within energy intervals (calculated as described in 
Section 3.3) corresponding to incident energies from 0.7 to 15 
MeV/nucIeon. 
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(4) For PHA events in the element bins from lithium through 

oxygen, a second count was obtained for the incident energy interval 

5. 9-9. 3 MeV/nucleon. The total energy boundaries for Li. Be and B 

7 9 11 

were ce’.culated using the nuclear masses of Li, Be and B 
respectively. 

The 0.7-15 MeV/nuclecn interval is optimum for the elements 
carbon throug.i nickel in the sen'^n that this is the three detector 
response interv '. common to these elements in all the LETs used. 
Likewise, 5. 9-9. 3 .MeV/nucleon is o^ mum for the elements lithium 
through oxygen. 

Three parameter analysis for iron required special attention as 
discussed in Section 3.5.2 and in Appendix B. 

3.5.1 The Charge Consistency Check 

Figure 3.7 shows a plot of Zl versus Z2 for a sample of events 
from LET C of Voyager 1. The events with consistent Zl and Z2 charge 
measurements fall in clumps along the diagonal, while other events 
fall in off-diagonal locations that are characteristic of various 
"background" effects to be discussed below. An expanded view of the 
events along the diagonal in Figure 3.7 is shown in Figure 3.8, a plot of 
(Z1-Z2) versus (Zl+Z2)/2. The lines indicate the consistency 
requirement used: 


Z1-Z2I < .0164(Z1+Z2) for (ZH-Z2)/2 > 5.5, 
< .180 for (Zl+Z2)/2 g 5.5 . 
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Figure 3.7 

A meatier plot of the charge measurements Zl and Z2 for PHA events 
obtained with LET C of Vo)'oger 1 during the period September 1977 to 
iMay 1970. (To prevent plot saturation, only every tenth event is 
plotted in the region where Zl and Z2 are both less than 9.) 




Figure 3.6 


A scatter plot of (Z1-Z2) versus (Zl+Z2)/2 for the same PHA events as 
in Figure 3.7. The lines indicate tiie consistency requirement used. 



2)/2 



In order to eliminate only a negligible fraction of "good" events this 
requirement was chosen to cut the distribution of (Z1-Z2) at about 
± 3 sigma for each element from carbon through nickel. The quality of 
the resulting data is shown in Figure 3.9, a (Zl+Z2)/2 histogram of 
the charge consistent events with energies from 8.7 to 15 MeV/nucleon 
acquired with the LETs on both Voyagers during the seven largest solar 
energetic particle events observed in the September 1977 to May 1978 
time period. The rms charge resolution ranges from 0.08 units at 
carbon to 0.27 units at iron. The charge boundaries used to sort the 
events into various element bins were chosen by inspection of Figure 
3.9 and are indicated by arrows, Due to spillover from the iron peak, 
abundances of Co and Mn were not obtained. Portions of the charge 
scale (near Co and Mn) that were not included in any other element 
bin are indicated by shading between the arrows of figure 3.9. 

The effect of the consistency check is shown by the comparison of 
Figure 3.9 with Figure 3.10, a Z1 histogram of all of the PHA events 
from the same time periods and telescopes as in Figure 3.9. Table 3.3 
compares the number of PHA events obtained for each element with 
and without the consistency check. This comparison indicates that, in 
the 8.7-15 MeV/nucleon energy interval, two and three parameter 
analysis yield similar (within 3 percent) abundances for the major 
elements (such as C, 0, Ne, Mg, and Si), but that three parameter 
analysis is necessary to reduce the background for the rarer elements 
(such as N, Na, AI, and Cr) which have more abundant neighboring 
elements of higher Z. 
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Figure 3.9 

A histogram of (Zl+Z2)/2 for charge consistent PHA events 
corresponding to nuclei with incident energies from 0.7 to 15 
MeV/nucleon. The histogram includes PHA events from all the LETs used 
on both Voyager spacecraft, summed over the seven major flare events 
in the September 1977 to May 1978 period. The data above oxygen and 
belc«v carbon are replotted with an expanded vertical scale. The arrows 
indicate the charge boundaries used to count PHA everts 
corresponding to the different elements. (The cross-hatched regions of 
the charge scale are discussed in the text.) 




Figure 3.10 


A histogram of the Zl charge measurement, including PHA events 
summed over the same time periods and LET telescopes as in Figure 
3.9. The effect of the charge consistency check may be seen by the 
comparison of this histogram, to the histogram of Figure 3.9. 
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TABLE 3.3 

The Effect of the Charge Consistency Check 


PHA Event Count 


z 

Element 

Charge Range 

without 

check 

with 

check 

Rejected 
Fracti on 

3 

Li 

2.5 


3.5 

838 

14 

0.983 ± 0.004 

4 

Be 

3.5 

- 

4.5 

413 

3 

0.993 ± 0.004 

5 

B 

4.5 

- 

5.5 

384 

2 

0.995 ± 0.004 

6 

C 

5.5 

- 

6.5 

8058 

7026 

0.128 ± 0.004 

7 

N 

6.5 

- 

7.5 

2729 

2096 

0.232 ± 0.008 

8 

0 

7.5 

- 

8.65 

19323 

17929 

0.072 ± 0.002 

9 

F 

8.65 

- 

9.5 

107 

4 

0.963 ± 0.018 

10 

Ne 

9.5 

- 

10.65 

2257 

2050 

0.092 ± 0.006 

11 

Na 

10.65 

- 

11.5 

308 

195 

0.367 ± 0.027 

12 

Mg 

11.5 

- 

12.6 

3564 

3299 

0.074 ± 0.004 

13 

A 1 

12.6 

- 

13.4 

328 

237 

0.277 ± 0.025 

14 

Si 

13.4 

- 

14.6 

3003 

2842 

0.054 ± 0.004 

15 

P 

14.6 

- 

15.4 

37 

10 

0.730 ± 0.073 

16 

S 

15.4 

- 

16.6 

613 

578 

0.057 ± 0.009 

17 

Cl 

16.6 

- 

17.5 

17 

10 

0.412 ± 0.119 

18 

Ar 

17.5 

- 

18.5 

68 

61 

0.103 ± 0.037 

19 

K 

18.5 

- 

19.4 

23 

13 

0.435 ± 0.103 

20 

Ca 

19.4 

- 

20.6 

201 

181 

0.100 ± 0.021 

21 

Sc 

20.6 

- 

21.5 

13 

3 

0.769 ± 0.117 

22 

Ti 

21.5 

- 

22.5 

25 

4 

0.840 ± 0.073 

23 

V 

22.5 

- 

23.5 

19 

3 

0.842 ± 0.084 

24 

Cr 

23.5 

- 

24.6 

80 

38 

0.525 ± 0.056 

26 

Fe 

25.25 

• 

27.0 

1879 

1570 

0.164 ± 0.009 

28 

Ni 

27.2 

- 

29.0 

115 

82 

0.287 ± 0.042 
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3.5.2 Dackground 

Most of the events eliminated by the consistency check are due to 
incomplete charge collection near the edge of the active areas of 
detectors Ll and L2. Edge effects in LI produce the events with 
deficient Zl charge measurements, but normal Z2 measurements, that 
are observed in Figure 3.7 in bands extending downward from the 
clumps along the diagonal. Edge effects in L2 result in the bands 
extending leftward from the clumps. (These bands are poorly defined 
and curve downward because both the Zl and Z2 charge measurement.s 
depend on the L2 pulse height.) The number of events in these bands 
indicates that the area of the "edge" of Ll or L2 is roughly independent 
of nuclear species and is about 6 percent of the fully active area. Edge 
effects are effectively removed by the consistency check since the 
probability for a heavy nucleus to masquerade as a lighter one by 
hitting the edge of both Ll and L2 is small. 

The band of events in Figure 3.7 which are near Z2=l are likely due 
to the accidental coincidence of a low energy heavy nucleus which stops 
in Ll and a proton which penetrates Ll and L2, and stops in L3. The 
smear of events near (Zl ,Z2 = 3) may be caused by either electronic 
pileup of protons and alphas or by nuclear interactions of these 
particles in the LET detectors or nearby material. 

In addition to the above background efl'ects there is an unusual 
one (to be discussed in detail in Appendix B) which resulted in the 
clump of events near (Zl=26 , Z2=28). These events have Ll and L3 
pulse heights which are appropriate for iron nuclei, but have L2 pulse 
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heights which are too large by about 10 percent. This effect caused a 
disproportionately large fraction of iron nuclei to be rejected by the 
consistency check (see Table 3.3). and required a correction of about 20 
percent to the 0.7-15 MeV/nucleon iron abundances as discussed in 
Appendix B. 


3.6 Three Parameter ”4-^ ”E) Analysis for Helium Nuclei 

dx dx 

The helium abundance and flux measurements presented in this 
work are based on three detector PHA events selected as follows; 


(1) Helium PHA events were selected from the LET Z< 3 events 
using a "slant" threshold on the L2 and L3 pulse heights: 




A sample of the selected events is shown in Figure 3.11(a), a scatter 
plot of L2 energy versus L3 energy. 

(2) For each event two mass measurements were calculated and 
required to be consistent (as described below). The rejected events of 
the sample shown in Figure 3.11a are replotted in panel (b), while the 
accepted events are shown in panel (c). Some of the rejected events 
lie on the L2 versus L3 helium track and are due to edge effects in 
LI. The number of these events indicates that for He nuclei the area 
of the edge of Ll is effectively 6 percent of the fully active area of Ll. 
This is the same percentage as found for the 3 nuclei, suggesting 

again that the edge effects are independent of Z and that their 
rejection by the consistency check does not introduce a 


Figure 3.11 


Scatter plots of the energy deposited in detector L2 versus the energy 
deposited in L3 for a sample of Z< 3 PHA events obtained with LET C of 
Voyager 1. Panel (a) shows candidate He PHA events selected using a 
"slant" threshold (indicated by the sloping line) to discriminate 
against proton PHA events (which fall below the line and are not 
shown). The candidate He PHA events which did not pass the mass 
consistency check (see text) are plotted in panel (b), while the 
accepted events are shown in panel (c). (Note that the scatter plots do 
not show the effects of discrete PHA channel numbers — a 
pseudo-random number with uniform distribution between 0.0 and 1.0 
was added to each pulse height channel number before conversion to 
energy units; see later text.) 
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Z-dependent bias. The rejected events which lie below the track may be 
caused by front incident protons which undergo nuclear interaction 
in L3 or by the accidental coincidence of a particle in LI with a side 
incident proton in L2 and L3. 

(3) Events with average mass measurement in the interval 3 to 5 
amu were sorted into various energy bins by their total energy 
measurement, Eli+Els+Eu. The energy bins were chosen to be 
within the response /ange of all the LETs used and correspond (as 
discussed in Section 3.4) to the incident energy/nucleon 
intervals listed in Table 4.4. 

The two mass measurements, Ml and M2, were calculated by 
solving: 


T'li “ R(El)+Ei^+El3 , Ml ) - R(Euj+E;j , M 1 ) 


and 


Tl2 = R(El2+E;j , M2) — R(El3 . M2) . 

Tli and Ti^ are the mean pathlengths for the appropriate Ll and L2 
detectors. Elj , Els ®^id El 3 are the energies deposited in detectors Ll, 
L2 and L3. R(E,M) is the range in silicon of a nucleus with nuclear 
charge Z = 2, and mass M, as scaled from the ^He range~energy relation 
Rhc of Vidor (1975) as follows: 

where Mhs is the mass of the “^He nucleus. 
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In order to properly account for the effects of finite PHA 
channel widths, a psuedo-random number with uniform distribution 
from 0.0 to 1.0 was added to each Ll, L2 and L3 pulse height channel 
number before the conversion to energies Ej^j , Elz and E^. Since the 
channel widths are small (about 70 KeV for detectors Ll and L2, 
and 500 KeV for L3) the approximation of a uniform distribution of 
energy measurements over each pulse height channel is accurate. 
The main effect of the random number addition is to produce 
continuous, rather than discrete, distributions for the mass 

measurments Ml and M2. 

The means of the Ml and M2 ^He mass distributions were found to 
vary from LET to LET, and as a function of total energy and time for a 
given LET, by less than about 0.2 amu. The means were typically within 
0.2 amu of the "^He nuclear mass (ft 4 amu) suggesting that 
systematic errors in the energy measurements Elj and Elz were less 
than about 2 percent for energies near 6 MeV, while systematic errors 
in the E^a energy measurements were less than about 3 percent near 30 
MeV. However, since 0.2 amu is a significant fraction of the helium mass 
resolution (ft 0.3 amu) achieved for individual LETs, the mass 
measurements Ml and M2 were adjusted (to obtain Mladj and M2adj) 
by the addition of correction factors (listed in Table 3.4) chosen for 
each LET to shift the mean masses to 4 amu. 

The mass consistency requirement: 

I ^ ladj ” I < 1.0 

was chosen to cut the (Ml^dj ~ M2adj) distributions at about ± (3 sigma 
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TABLE 3.4 

LET Helium Mass Correction Factors 


Voyager 1 

LET A 
LET B 
LET C 
LET D 

Voyager 2 

LET A 
LET B 
LET C 
LET D 


ami 


AM2 


0.216 

0.211 

0.197 

0.138 

0.123 

0.110 

0.047 

0.004 

0.209 

-0.050 

0.102 

-0.058 

-0.007 

-0.045 
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+ .25 omu) In order to reject only a negligible fraction of good events, 
The mass resolution achieved is Illustrated in Figure 3.12 which shows 
two (Ml,,jj+M2,dj)/2 histograms of mass consistent PHA events from 
the LETs on both Voyagers. Histogram (a) is for the flare event of 

r| 

September 9, 1977, while histogram (b) is for the small He rich event of 
October 13, 1977. Since no He peak was seen in the mass histograms 
for any of the large flare events studied in this paper, the mass 

selection interval 3 to 5 amu was chosen to include essentially all good 

4 3 

He events, with negligible He contamination. 

In one of the major flare events studied (the April 29, 1978 event) 
the He mass resolution was degraded due to high count rates. In this 
case, the consistency check was relaxed to i M - M2ujj | < 2.0 and 
the mass selection interval was increased to 2-5 amu to prevent the 
rejection of "^He PHA events which, were apparently affected by baseline 
shifts in the LI pulse height analyzers. 


3.7 Fluence and Relative Abundance Measurements 
3.7.1 Definitions and Terminology 

dJz 

The 'differential flux" of nuclei with nuclear charge Z is 

dE 

defined, as a function of energy/nucleon E, direction fi , location r. and 
time t, by: 


d J, 

--- (E , n , r , t) 5 N/ (dA • dfl • dt • dE) 
dL 


particles/ (cm^ sr sec MeV/nucleon) 
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Figure 3.12 

Helium mass histograms, summed over mass consistent PHA events 
from the LETs on both Voyager spacecraft. Histogram (a) is for the 
flare event of September 9, 1977; (b) is for a small *^He~rich event of 
October 13. 1977. 


COUNTS PER BIN 
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where, in the time interval dt (sec) containing t, N is the number of 
nuclei which (a) cross a surface area dA (cm®) that is oriented normal to 
Q , (b) have trajectories within the solid angle dfi (sr) centered on f) , and 
(c) have kinetic energies in the range dE about E (MeV/nucleon). 

The 'Integral flux" Jj is defined as the integral with respect to 
energy/nucleon of the differential flux, and therefore refers to the flux 
contained within some specified °nergy/nucleon interval. 

The integration with respect to time of a flux yields a "fluence". 

The SEP "relative abundance" of two different nuclear species is 
defined here as the ratio of their differential fluxes at common 
values of E, n . r and t. 

Generally, measurements of SEP flux, fluence, and relative 
abundances are averaged over some range of energies, incidence angles 
and times and are performed at a particular location in the 
heliosphere. The averaging techniques employed in the fluence and 
relative abundance measurements presented in this thesis are described 
below. 

3.7.2 Fluence Measurements 

For each spacecraft, the integral flux of nuclei, with nuclear 

charge Z, in energy interval i, was calculated for each three-hour 
interval j by averaging the fluxes measured in each LET k as follows: 






m 
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u " NPHAj,u ' " TiM 
particles/ (cm^ • sr • sec) 


(3.2) 


The terms are explained below: 

• K is the number (either 3 or 4) of LETs used; LET B on 
Voyager 1 and, at times, LET C on Voyager 2 were excluded due 
to the problems discussed in Appendix B. 

• ^2 j j j, is the number of PHA events (selected by two or three 
parameter analysis as discussed in previous sections) for 
element Z , energy interval i , three-hour interval j , and LET 
k . 

• NPHAjk is the total number of PHA events read out from the 
appropriate (Z< 3 or Z^ 3) event buffer for LET k during the 
three-hour interval j . 

• Cj is the number of Z< 3 or Z^ 3 counts recorded by the rate 
accumulator for LET k during interval j . 

• Tj It is the length of time (in seconds) during which the 
counts Cj It were accumulated. (This time was often less than 
three hours sine the spacecraft were not continuously 
monitored.) 

• LTj.u is a small correction for dead time of the rate 
accumulator. Laboratory measurements (Gehrels, private 
communication, 1980) of the high count rate response of the 
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LET electronics indicate that: 

LT,.k = (l -T1 •Rlj.kKl -T2 R2n) 

where Rlj,k and R4j_j, are the count rates for detectors LI 
and L4, with effective dead times per count of Tl = 18 /isec and 
T4 = 20 fxacc respectively. 

• is the geometry factor (in cm^ steradian) for LET k (from 
Gehrels 1980). 

For each spacecraft, the integral fluence measurement , 

summed over a flare period t , including three-hour intervals indexed 
by the subscript j , was calculated using; 

1 1 E Jz i j ■ ^T particles/ (cm^ • sr) 

J 

where AT = 10000 seconas (3 hours) and the subscript i still indicates 
the energy/nucleon interval. Since the spacecraft were not continuously 
monitored, some three— hour intervals were necessarily excluded from 
the sum. For each flare event a rough estimate (included in Tables 
4.2 through 4.7) of the missing fluence was made by interpolating the 
counting rate for ZS 3 nuclei through data gaps. 

The statistical uncertainty in the flux and fluence 
measurements is dominated by the uncertainty in the number of PHA 
events (Nzjj ,k) obtained for a particular element and energy bin. N 2 ,ij,i< is 
always much smaller than either the total number of PHA events 
(NPHAj.,) or the number of rate rcalar counts (Cji<). Thus, the 
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statistical uncertainty was estimated using the right hand side of 

j. 

equation (3.2) with Nz j jjt replaced by ® . 

3.7.3 Relative Abundance Measurements 

Two different types of relative abundance measurements were 
made for each flare event: "flux weighted" and THA event weighted". 
The flux weighted relative abundance measurement (A 2 i,z 2 ,i^F,t 
nuclear species Zl and Z2, for the common energy/ nucleon 
interval i. and time period t, is defined as the ratio of fluences; 

^Zl,Z2.1^F,t = " zi.i.t /^'^Z2.i,t • 

The PHA event weighted relative abundance measurement Azi_z2.i P.i 
may be defined for nuclear species Zl and Z2 which are pulse height 
analyzed with the same priority (i.e. Zl,Z2< 3 or Z1-Z2 ^ 3). as the ratio 
of PHA event counts: 

^Azi,z2,i)p.t = Nzu.t /Nz2,i.t • 

If the true relative abundance of the two nuclear species is constant 
during the averaging time, t, then (A)p and (A)p will be the same within 
statistical error. However, if the true relative abundance varies 
during t and if the flux is high enough to saturate the telemetry rate 
for PHA events, then the two measurements may be systematically 
different. In the next chapter the time variation of relative abundances 
during individual flare events will be examined, and "flux" and 'PHA 
event" weighted measurements will be compared. 
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It is important to note that for a flare event in which the PHA 
event telemetry saturates, the flux weighted abundance 

measurements will generally have larger statistical uncertainties than 
the corresponding PHA event weighted measurements. This is because a 
flux weighted measurement may depend primarily on a 

disproportionately small number of PHA events acquired during the 
time of peak flux, while 'PHA event weighting" makes full use of PHA 
events acquired throughout the averaging time period, giving each 
event equal statistical weight. Since PHA event weighting produces 
smaller statistical uncertainty, it was desirable to extend the technique 
to obtain relative abundances of nuclei which are not pulse height 
analyzed with equal probability. In particular, the PHA event 
weighted abundance measurement A^e.z.! Pt helium relative to 
element Z 3) for a common energy/nucleon interval i, and a time 
period t, was defined separately for each spacecraft by : 


^Ahc.z.i^p.i = I 06 • exp 





where: 


• and Nz.ij are the numbers of PHA events for helium 
and element Z acquired during a time subinterval j (of length 
chosen to be six hours) and summed over the LETs used. 

• PnejandPzj are the pulse height analysis probabilities for 
Z< 3 and ZS 3 nuclei respectively 


PHe,j ^ NPHA(2<3 )j /C(ZO) j , 
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Pzj * NPHA(25 3 )j / C(ZJ 3)J . 

where NPHA and C are the numbers of PHA events and rate 
scalar counts (for Z< 3 or 3 nuclei) obtained during time j, 
and summed over the LETs used. 


• wj is a weighting factor chosen to minimize an estimate Cp of 
the statistical uncertainty in Ahc.z.i P.t • 




He.i.j 


+ 


1 



(Sw,) 

J 




P.t 


• The factor of 1.06 accounts for the 6 percent difference 
between the effective geometry factor for He nuclei (which 


d E d E 

were counted by three parameter E analvsis) and 

dx dx 

the effective geometry factor for the higher Z nuclei (which 
for this measurement where counted using two parameter 


dx 


-E 


analysis). 


• The logarithmic average was used to obtain the desirable 
property, (Ahc,z>P = V- 
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Chapter 4 
Observationa 


4.1 Overview 

The seven largest solar energetic ’ rticle (SEP) events in the 
September 1977 to May 1978 period are selected for study. The time 
development of each event is examined in order to select an optimum 
averaging time period for SEP relative abundance measurements. 
Abundance measurements of nuclei from He (Z=2) through Ni (Z = 28) 
are presented for each flare event and are compared to other 
reported measurements which are available for four of the seven flare 
events. Fluence measurements are presented for the more abundant 
elements (He, C. 0, Ne. Mg, Si, and Fe + Ni) in several energy bins from 3 
to 30 MeV/nucleon. The possible systematic errors in the presented 
measurements are estimated. 

' 7 Energetic Particle Event Selection 

Figure 4.1 shows the hourly averaged flux of Z> 3 nuclei with 
energies above about 3 MeV'/nucleon, as measured by the LETs on 
Voyager 1 during the period September 1977 through May 1978. The 
seven largest energetic particle events (labeled i through vii) were 
selected for study because of the statistical accuracy available for Z§ 3 
abundance measurements. Optical, radio and/or X-ray data (CofTey 
1977, 1978) indicate that fnr each of the first four events a large solar 
flare occurred within a few hours before the arrival of relativistic 
electrons (as measured by the Voyager High Energy Telescopes). The last 
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Figure 4.1 

The hourly averaged flux of 3 nuclei measured with the LETs on 
Voyager 1. The energy threshold for detection varies from about 3 
MeV/nucleon for carbon nuclei, to 5 MeV/nucleon for iron nuclei. The 
seven largest energetic particle events are labeled i through vii. (The 
quantization of the flux measurement at low levels corresponds to the 
detection of a small integer number of nuclei.) 
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three events were not well separated, but were part of a general intensity 
increase accompanied by several large solar flares. 

The selected energetic particle events are large: The peak flux of 
protons (>20 MeV) observed in the events ranges from about 20 to 500 
particles/(cm* • sr • sec) placing all seven events in the upper 5 percent 
of the size distribution for energetic particle events recorded near the 
time of the previous solar maximum. 1967 through 1969 (Svestka 
1976). The first three of the seven events were observed by neutron 
monitors at ground level (CofTey 1977, 1970). 

A compilation of references to some recent energetic particle 
observations of the first four events may be found in Wibberenz 
(1979). The importance, location on the sun, and time of maximum for 
the large optical flares preceding energetic particle events i, iii and iv 
are respectively: 3B, N08W57, 1977 September 19 10:30 UT: 2B, 

N24W40, 1977 November 22 10:05 UT; and 2N. N 1 5W20, 1 9?0 February 
13 1:43 UT (Co^Tey 1977, 1976) Energetic particle event ii probably 

originated in a flare which occurred just beyond the west limb in the 
same active region which produced the flare of 1977 September 19. 
(Type II and III radio emissions were observed at about 5:55 UT of 1977 
September 24; CofTey 1976) 

4.3 Energetic Particle E .’ent Time Profiles 

The time structure of the seven selected energetic particle events 
is displayed in Figures 4.2 through 4.6. Shown separately for each 
Voyager spacecraft are: (1) counting rates corresponding to protons 


Figures 4.2 through 4.6 


Energetic particle event time profiles. Various counting rates 
(labeled LAI, BSP, and ZZ 3), element ratios (C/0, Ne/0, Mg/0, Si/0 and 
[Fe+Ni]/0), and the oxygen spectral index (7oxygen) ®^e averaged in 
successive three hour intervals and are plotted separately for 

Voyager 1 ( — ,•) and Voyager 2 (• ,o). The counting rates LAI, BSP, and 

ZZ 3 correspo.;-! respectively to protons above about 0.5 MeV, protons 
above 20 MeV, and ZZ 3 nuclei above an energy threshold which varies 
from about 3 MeV/nucleon for carbon to 5 MeV/nucleon for iron. The 
element ratios are for the common interval 5-15 MeV/nucleon, except 
iii the case of the He/0 ratio where the energy range is 4. 0-7. 8 
MeV/nucleon, The oxygen spectral index was computed from the ratio 
of oxygen PHA events obtained in the 4.0-6. 1 and 6.1-15 MeV/nucleon 
energy intervals. (The cross-hatched boxes are discussed in the text.) 
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above 0.5 MeV. protons above 20 MeV, and 3 nuclei above about 3 
MeV/nucleon, (2) the abundance ratios C/0, Nc/0, Mg/0, Si/0 and 
(Fe+Ni)/0 measured over the energy range 5.0-15.0 MeV/nucleon, (3) 
the He/0 ratio measured for the 4.0-7.f' MeV/nucleon energy 
interval, and (4) the spectral index y (assuming dJ/dE a E"’’) of 
oxygen, obtained from the ratio of oxygen PHA events in the 
6.1-15.0 and 4.0-6. 1 MeV/nucleon energy intervals. The main 
features of these plots are summarized below. 

(1) Generally, there is good agreement between the Voyager 1 and 
2 measurements of abundance ratios and the oxygen spectral index. 
Thus, for the relative abundance results presented later we combined 
the data from Voyagers 1 and 2 (improving statistical accuracy and 
averaging over any small spatial inhomogeneities which may exist 
between the spacecraft), 

(2) The (Fe+Ni)/0 ratio decreased by a factor of three to five during 
the first 12 to 24 hours of each of the energetic particle events which 
showed velocity dispersion. Velocity dispersion was identified in the 
first four events (Figures 4.2 through 4.4) by the earlier arrival of 
protons above 20 MeV relative to protons above 0.5 MeV and by a 
softening of the oxygen spectra early in the events. Time variations 
associated with velocity dispersion were also seen for the He/0 ratio, 
which increased during the 1977 November 22 and 1970 February 14 


events. 
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(3) In contrast to { 2 ), systematic time variations were not 
apparent for the C/0, Ne/0, Mg/0 and Si/0 ratios, which were constant 
within about ± 30 percent throughout each of the energetic particle 
events. 

(4) The last three events showed more complex time structure 
than the first four well-separated events. It is interesting to note that 
the 197B April 29 event showed no velocity dispersion, had a 
symmetrical intensity rise and fall, and occurred at Voyager 2 six 
hours earlier than at Voyager 1. Since Voyager 2 was located about 
1/0 AU sunward of Voyager 1, these features suggest that the energetic 
particles were trapped in a structure which was moving radially away 
from the sun at about 000-1000 km/sec. 

Abundance time variations, such as in (2) above, which are 
associated with velocity dispersion have been seen at somewhat lower 
energies by a number of experimenters (see, e.g.. Van Allen, 

Venketarangan and Venkatesan 1974; O'Gallagher et. al. 1976; 
Scholar et. al. 1970; von Rosenvinge and Reames 1979) and are thought 
to be a propagation effect in which nuclei with the same velocity, but 
different charge to mass ratios, diffuse from the sun to the 

spacecraft at different rates. O'Gallagher et. al. (1976) studied the 
arrival time distributions of H, He, C, 0, and Fe nuclei (0,5 to 5 

MeV/nucleon) in one flare event which showed particularly regular 
intensity versus lime profiles and concluded that the H through 0 nuclei 
were almost fully stripped but that the Fe nuclei had retained about 
16±5 of their 26 electrons (charge state of 10 ± 5), Although 

direct charge state measurements at these energies are not available 


for comparison, these results are consistent with measurements at 
lower energies by Sciambi et. al. 1977 and Gloeckler et. al. 1976 who 
found mean charge states of 5.8 for carbon, 6.3 for oxygen (averaging 
over nine flares, for energies 15-600 KeV/nucleon), and 11.6 for iron 
(from a single event, for energies 6-250 Kev/nucleon). 

The significance of the observed abundance time variations will 
be further discussed in the next chapter. Here we consider the effect of 
the time variations on time averaged abundance measurements for 
individual flare events. Abundance time variations have been studied 
recently by Scholer et. al. (1976) using a relatively standard 
propagation model which incorporates rigidity dependent diffusion, 
convection and adiabatic deceleration. The results are in qualitative 
agreement with the observed time dependence of the (Fe+Ni)/0 ratio, 
and suggest that the ratio observed later in an event, after the initial 
decrease of the (Fe+Ni)/0 ratio, approximates the injection ratio at the 
sun. This was taken into account in the selection of averaging time 
intervals (indicated by the shaded boxes of Figures 4.2 through 4.6) for 
the abundance measurements presented in the next section. 

For the four flare events which showed velocity dispersion (events 
i, ii, iii, and iv) the selected averaging periods exclude times early in 
the events when the Fe/0 ratio was decreasing. For the other three 
events which show more complex time structure the averaging time 
periods were selected differently. The averaging time for event v was 
chosen to include the period late in the event when the Fe/0 ratio was 
roughly constant. The averaging time tor event vi includes essentially 
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the whole event. For event vii, the avereging period was selected to 
avoid contamination by event vi. 

The sensitivity of the abundance measurements to the choice of 
averaging time period was checked by comparing Fe/0 
measurements averaged over the selected time periods to the 
corresponding measurements averaged over entire flare events. For 
those four flare events with the largest Fe/0 time variations, the results 
are listed in Table 4.1. In addition, the two averaging techniques "flux" 
and "PHA event" weighting (discussed in Section 3.7) are compared. For 
the Fe/0 relative abundance measurement, the effects of averaging 
technique and time period selection are seen to be comparable in size 
to the statistical errors, but are small compared with the flare to 
flare differences. 

4.4 Relative Abundances 

Abundance results for each of the seven energetic particle events 
are presented in Tables 4.2 and 4.3. The main results are the relative 
abundances of the elements carbon through nickel which are based 
entirely on three detector PHA events (see Section 3.5 and Appendix 
B) and were obtained for a single common velocity interval 
corresponding to the interval 8.7-15.0 MeV/nucleon. The He/Si 
ratios are for the LET response interval 4. 6-7.8 MeV/nucleon which is 
common to He and Si, and were obtained using three detector helium 
PHA events (see Section 3.6) and two detector silicon PHA events (see 
Section 3.4). All of the abundances are 'PHA event" averaged (see 
Section 3.7.2) over time periods chosen as discussed above, 
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Table 4.1 

Fe/0 Abundance Ratio 
(8.7-15 MeV/Nucleon) 




Solar Energetic 

Particle Event 



1977 

1977 

1977 

1978 


Sep. 19 

Sep. 24 

Nov. 22 

Feb 

. 13 

Full Event 






Flux Weighted 

0.28-0.05 

0.78±0.06 

0.31±0.04 

0.12 

-0.01 

PHA Event Weighted 

0.33:0.04 

0.78-0.05 

0.35±0.04 

0.13 

-0.01 

Partial Event 






(onset excluded) 






Flux Weighted 

0.24-0.04 

0.67±0.07 

0.26±0.05 

0.10 

-0.0] 

PHA Event Weighted 

0.29-0.04 

0.67-0.06 

0.30±0.05 

0.11 

±0.01 
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For most elements the abundances are normalized to silicon to 
facilitate comparison to solar abundances measured by other means. 
However, the elements Li, Be, and B have very low abundances in the sun 

4 * 

and their presence in our data would indicate spallation of energetic C, 
N and 0 nuclei. Thus, the Li, Be and B^*abundance upper limits are 
referenced to 0 and are based entirely on three parameter analysis in 
the LET response interval 5. 4-9, 3 MeV'/nucleon which is common to 
the elements Li through 0 (see Section 3.5), 

The listed (± lo) uncertainties include only the effect of 
counting statistics, which is generally large compared to the possible 
systematic errors (discussed below). Eighty-four percent confidence 
upper limits are quoted if no counts were obtained for an element or if 
the element showed no clear peak in the charge histogram, Figure 3,9. 

4.4.1 Systematic Error 

The main sources of systematic error in the relative 

abundance measurements of C through Ni (8.7- 15.0 MeV/nncIeon; Table 
4.2) arc: 

(l) Energy threshold errors. The actual energy threshold may 
differ slightly from 8.7 MeV/nucleon, and more importantly, may vary 
with nuclear charge Z as a result of (a) possible systematic errors in the 
energy calibration of the LET detectors and/or (b) deviations of the 
thickness of the thin A1 windows from the 3 ^m value specified by 
the manufacturer. The variation of the threshold with Z would arise 
since, near the 8.7 MeV/nucleon threshold, nuclei with different Z 
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deposit their energy In different proportions among the A! entrance 
window and the detectors LI. L2 and Lb. 

The size of the possible systematic error in the energy 
calibration of the LET detectors was checked in two different ways 
during the data analysis and was found to be about 2 percent for the Ll 
and L2 detectors (see Section 3,6) and about 3 percent for the L3 
detectors (see Sections 3.3 and 3.6). As a result of the manner in 
which energy is distributed among the detectors by incident nuclei 
near the 0.V MeV/nucleon threshold, a possible variation of the 
threshold with Z is most sensitive to a relative calibration error 
between the L3 energy and the average of Ll and L2 energies. This 
relative error was conservatively taken to be 5 percent, and the 
resulting Z-dependent shift of the 8.7 MeV/nucleon energy threshold 
was calculated for each element, carbon through nickel. Assuming 
typical energy spectra, dJ/dE a E"^ (see Section 4.5), the energy 
threshold shifts yield an error of about 0.5- |Zl-Z2i percent in the 
relative abundance of elements with nuclear charge Z1 and Z2, For 
adjacent elements this error is very small, and is only about 10 percent 
for the most widely separated elements carbon and nickel. Such 
possible errors are much smaller than the observed hare-to-flare 
variation (e.g. the Ni/C ratio measurements range over a factor of ten; 
see Table 4.2) and are usually smaller than the statistical error. 

The effect of possible A1 window thickness errors is also small. 
For example, in order to induce a Z-dependence La the relative 
abundance measurements of about 0.5- !Z1-Z2i (again, for energy 
spectra dJ/dE o E“^) a window would have to be twice as thick as the 


specitied 3 fim. 


(2) Charge assignment error and the '\ron" problem for Cr and Ni. 
For the elements with clear charge peaks in the charge histogram, 
Figure 3.9, the biases introduced by the data selection and charge 
assignment techniques are negligible (see Section 3.5), with the 
possible exception of the iron group nuclei. The background effect 
discussed in Appendix B necessitated corrections of about 20 percent 
to the iron abundances. The same background effect probably also 
occurs for Cr and Ni nuclei, however corrections were not applied for 
these elements. The error thus induced (about percent) is probably 
comparable in size, but opposite in sign, to that produced by 
spillover of Fe nuclei into the Cr and Ni peaks. Therefore a rough 
estimate of the systematic uncertainty in the the abundances of Cr and 
Ni relative to Fe is ± 20 percent. 

The main systematic error in the He/Si ratios (4. 6-7. 8 
MeV/nucleon; Table 4,3) results from the possible energy threshold 
errors discussed above. A 5 percent relative error in energy 
calibration between detectors L3 and L1+L2, together with E"^ energy 
spectra, gives an error of 11 percent in the He/Si ratio. The He/Si 
measurement in event vi may have an additional systematic error (on 
the order of 20 percent) due to the high count rate degradation of He 
mass resolution mentioned in Section 3.6. Note that this possible 
additional error is still small compared to the factor of three dip in the 
He/Si ratio during event vi (see Figure 4.6). Since both the possible 
systematic error and the time variation of the He/Si ratio in event vi are 
much larger than the statistical error, the statistical error is not 
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quoted in Table 4.3 and the reader is referred to this discussion. 

4.4.2 Comparison to Other Reported Observations 

The abundance results of Tables 4.2 and 4.3 have been compared to 
other reported observations — McGuire, von Roserivinge and 
McDonald (1970) for flare events i, ii, iii, and iv, and Dietrich and 
Simpson (197b) for event ii. For most of the elements (namely He, C. N. 
0, Ne, Na. Mg, Al, Si. S, Ar, Ca, and Fe) the various measurements 
generally agree within statistical errors. The largest 
non-statistical difference occurs for the Fe/0 ratio in event ii — we 
obtained Fe/0 = 0.70± 0.06, while McGuire, von Rosenvinge and 
McDonald (1979) found (Fe+Ni)/0 = 1.17± 0.11. This difference may 
be due to McGuire et. al. including the event onset period in their 
averaging time. 

However, our results do not support the high abundances of the 
rare elements B, F and Cr reported by Dietrich and Simpson (1978) for 
event ii. Their finite results are significantly higher than our upper 
limits for B and F and our finite value for Cr (as reported earlier in 
Cook et. al. 1979). The discrepancies may result from the fact that 
Dietrich and Simpson’s results were obtained at energies above 25 
MeV/nucleon, v;hile our data are for the lower energy interval 8.7-15 
MeV/nucleon. However, the B, F and Cr results of Dietrich and Simpson 
(1978) are based on only a few PHA events, so there is ccn'..e."n 
about contamination from the more abundant elements C, 0, Ne and Fe. 
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Dietrich and Simpson (1978) compared their SEP abundances to 
"solar system" abundances which were a mixture of results from 
different sources and concluded that (a) the enhancement of even-Z 
nuclei increased monotonically with Z and that (b) the rare nuclei (I:), F, 
Na, A1 and Cr) were additionally enhanced due to fragmentation of 
heavier nuclei in traversing *«0.6g/cm^ of solar atmospheric material 
(see Figure 1.3). We disagree with both conclusions. The 
over-enhance.ment of Na and Al. as well as the monotonicity of the 
enhancement of even-Z nuclei disappears if the SEP abundances are 
compared to a consistent set of abundances based only on solar- 
spectroscopy (see Cook et. al. 1979, McGuire, von Rosenvinge and 
McDonald 1979, and Chapter 5 of this thesis). For event ii. our 
abundances for Na, Al and Cr and upper limits for B and F are 
consistent with negligible matter traversal and with the upper limit of 
0.06 gm/cm^ reported by McGuire, von Rosenvinge and McDonald (1979). 

4.5 Flucnce Measurements and Energy Spectra 

Fluence measurements of He, C, 0, Ne, Mg, Si and Fe nuclei for 
several energy/nucleon bins are presented in Tables 4.4 through 4. 10; 
one table for each energetic particle event. The fluence measurements 
were integrated (as described in Section 3.7) over time periods chosen 
to include e»jtire energetic particle events, and are averaged over 
Voyagers 1 and 2. Table 4.11 lists the fluences measured for the entire 
time period from September 1977 to May 1978. 


TABLE 4.4 

Fluence Measurements for the Flare Period: 1977 September 19 1200 UT to 1977 Septeater 22 1500 UT 
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Fluence Heasurpn^e^its for the Flare Period: 1977 Novetnber 22 1200 UT to 1977 Novert)er 27 0300 UT 
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Huence Measurements for the Entire Period: 1977 SepteiMier 19 1200 UT to 1978 Hay 7 1800 UT 


112 





1/1 re 01 

U7UJ L. >K U 

H o c 

80 o CT ^ 


. O t) 

W O) 

• o 0) T) 

> 0) O 

uj ui o cr 3 

• pv re CT 

• ■ W 01 < 

6 > X ^ 

E re 4-> 01 

4/1 01 W U 

40 — C O -e- 

o o 01 re c 

bj uj re 1/1 c t/i 

lTt«3 01 0i 0; 

^ 1/9 u u T? 

• • ^ C U 3 

*o »> C o 07 »— 

2 3 a. u 

o e- c 

o Lu ^ 



113 


The fluence measurements presented here are intended to serve 
several purposes. In Chapter 5 they will be used to examine the 
energy /nucleon dependence of SEP composition in the 5 to 15 
MeV/nucleon interval. In addition, the fluence measurements should 
eventually be combined with forthconning results from the Voyager Low 
Energy Charged Particle Experiment and the Voyager High Energy 
Telescopes to study SEP composition over a wide energy range from 
below 1 to above 100 MeV/nucleon. Finally, the fluence measurements 
presented here are important to a number of astrophysical problems 
which are not addressed in this thesis, including the study of the 
effects of SEPs on lunar rocks and soil (see, e.g.. Price et. al. 1974) 
and the consideration of flaring stars as possible injectors of energetic 
particles into the galactic cosmic ray accelerator (see. e.g.. Montmerle 
1979). 


In anticipation of the discussion in Chapter 5, the fluence 
measurements of the elements C, 0, Ne, Mg. Si. and Fe-fNi in each 
flare event were fit by power law functions of the form A ■ (E/'7.5)"’' by 
minimizing the function: 

xHa.7) = S '"a • (E/7.5)-^dEjVa,2 (4.1) 

where Wj and are the fluence and its statistical error for 

energy/nucleon bin i having lower and upper energy/nucleon 
boundaries Ej and E( 4 .j respectively, The same six energy/nucleon bins 
in the interval 5-10 .MeV/nucleon were used for each element. The 
fluence measurements for these energy/nucleon bins and the best fit 
power law functions are plotted in Figure 4.7, The best fit spectral 


Figure 4.7 


Energy spectra of the elements C, 0, Ne, Mg, Si, and Fe-t-Ni for the flare 
events i through vii. Also shown are best fit power law functions of 
the form A(E/7.5)”’', where E is kinetic energy per nucleon, (x® is 
defined in the text.) 
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indices are suimnarized in Table 4.12. 


The normal.'zing energy, 7.5 MeV/nucleon, was chosen such that A 
and y were roughly independent fitting parameters, allowing the 
statistical uncertainty in A to be estimated by: 





1 

2 


The statistical uncertainty in y was estimated (as described by 
Bevington 1969, p. 243) as that increment to the best fit value of y which 
increased ^ by 0.25 fr im its minimum value. 


Note that the above spectral indices are derived from fluence 
measurements which are integrated over entire SEP events, rather than 
over the time periods selected in Section 4.3 to exclude event onset 
times. Due to the softening of energy spectra which occurs throughout 
some events (see, e g., the oxygen spectral index in Figure 4.4) the 
spectral index measurement for a given element typically does depend 
somewhat on the choice of averaging time period. However, in the 
discussion in Chapter 5 we will be interested only in the differences 
between the energy spectra of the different elements. These 
differences were found to be insensitive to whether or not the event 
onset times were included in the averaging time periods. 


4.5.1 Systematic Error 

The relative systematic errors between fluence measurements for 
the different elements and energy bins are dominated by LET 
detector energy calibration errors (except in the case of He in event vi, 
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I 
r 

I 

as Indicated later in this section and in Section 4.4.1) The effect of 
such possible errors on relative abundance measurements (in 
I particular, for the 0.7-15 MeV/nucleon interval) was discussed in 

section 4.4.1, and found to be small; about 0.5 ■ |Z1-Z2| percent for the 

I relative abundance of elements with nuclear charge Z1 and Z2. Here we 

I 

I consider the effect of possible energy calibration errors on the spectral 

i 

I index 7. The possible error in 7 was estimated by (1) calculating the 

i shift of energy thresholds (Ej, see above section) induced by a 5 percent 

relative error between the L3 and L1+L2 energy measurements and (2) 

I using the shifted Ej in expression 4.1 to recompute 7. The average 

* percentage changes in 7 for the elements C, 0, Ne, Mg, Si, and Fe were 

about 7, 4, 1.6, 1.2, 0.5 and -1.0 respectively. While these possible 
errors are small re-ative to many of the observed differences in spectral 
index among the different elements (see Table 4.12), they are sometimes 
comparable to the statistical errors and therefore are taken into 
account in the spectral index comparisons of the next chapter. 

The main sources of possible absolute error in the fluence 
measurements of Tables 4.4 through 4.1 1 are discussed below; 

(1) Geometry factor. The uncertainty in the absolute LET 
geometry factors is estimated as 6 percent since the area of the "edge” 
of detectors Ll and L2 was found to be about 6 percent of their 
active areas (see Sections 3.5.2 and 3.6). 


(2) Energy calibration. Conservatively estimating the absolute 
error of particle energy measurement at 5 percent and assuming a 
typically steep energy spectrum (dJ/dE a E"^) gives a fluence error of 15 
percent. 

(3) Data gaps. The fluence missing due to data gaps was 
estimated for each flare event and ranges from 0.5 to 35 percent of the 
measured fluence. Thus, the error in the uncorrected fluence 
measurements (i.e. those listed in Tables 4.4 through 4.11) due to data 
gaps is S 35 percent. Since the gaps were generally short (about 3 to 12 
hours) and the flux was typically a fairly smooth function of time, the 
estimates of missing fluence are probably accurate to better than about 
± 50 percent. Thus, if the fluence measurements were corrected for the 
missing fluence, the residual uncertainty due to data gaps would range 
from about 0 to 20 percent for the different flare events. 

(4) Anisotropy. The fluence measurements of Tables 4.4-4.11 may 
differ from omni-directional fluence measurements as the result of 
particle anisotropy. For each flare event, the size of such possible 
differences was estimated by comparing (a) fluence measurements 
averaged over all four Voyager 1 LETs to (b) fluence measurements 
averaged over only Voyager 1 LETs A and C, which, pointing in 
opposite directions, measure the omni-directional component of the 
fluence. The differences between the fluence measurements (a) and 
(b) ranged from about 0 to 20 percent with a typical value of about 10 
percent. 


(5) High count rates. Possible systematic error due to high 
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counting rates only occurs for He in event vi, where as mentioned 
earlier, the mass resolution was somewhat degraded, Although for 
event vi, the He moss comsislency check was relaxed and dead time 
corrections were applied, the residual systematic error in the He 
fluence measurements may be on the order of 20 percent. 

Combining the above error estimates gives an absolute error for 
the fluence measurements of Tables 4,4 through 4,11 in the range of 
about 15 to 50 percent, depending mainly on the data gaps. 


121 


Chapter 5 
Discussion 


5.1 Overview 

With the goal of understanding the relationship between SEP 
elemental composition and the composition of the sun, this chapter 
discusses the systematics of the SEP observations presented in 

Chapter 4. As discussed in the Introduction, we will (a) examine 
the time and velocity (energy/nucleon) dependence of the SEP 
composition in the individual flare events, (b) select for further 
study those flare events which have SEP composition that is 

approximately independent of energy/nucleon, (c) compare the 
selected SEP composition results from flare to flare, and (d) compare 
the selected SEP composition results to solar abundance 
measurements from other sources: spectroscopy of the photosphere 
and corona, and solar wind measurements. 

An overview of our SEP abundance results and their 
relationship to other solar abundance measurements is shown in Figure 
5.1. a plot of the SEP abundance measurements (from Tables 4.2 and 

4.3) for all seven flare events, along with abundances for the 

photosphere, corona and solar wind. Although the SEP composition is 
seen to vary from flare to flare, the average SEP abundances, when 
normalised to silicon, are similar to abundances from the other 
sources for all the elements shown except C, N’, and 0, where the SEP 
values are persistently lov.- relative to the photosphere. For oxygen, the 
SEP abundances are similar to the solar wind and active coronal region 


Figure 5.1 


Elemental abundances relative to silicon. References: ( ) Meyer and 
Reeves (1977), (t ) Same et. al. (1975), { § ) Parkinson et. al. (1977). 
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values. 

In the following discussion of the time and velocity 
dependence of SEP abundances in individual flare events, we are 
particularly interested in whether or not there is an indication of SEP 
acceleration and/or propagation bias which may account for the 
persistent differences seen between SEP elemental composition and the 
results of photospheric spectroscopy. 

5.2 Time and Velocity Dependence of SEP Abttndances for Individual 
Flare Events 

The time development of the seven SEP events was discussed (in 
Section 4.3) in connection with the selection of averaging time 
periods for the SEP abundance measurements fo ' each flare event. 
Clear evidence for systematic abundance time variations was seen for 
the Fe/0 ratio (5-15 MeV/nucleon), which decreased by a factor of 
three to five in each of the four flare events that showed velocity 
dispersion in the particle arrival times. In two of these four events the 
He/0 ratio (4. 0-7. 8 MeV/nucleon) increased with time. However, the key 
point for the discussion in this chapter is that systematic time 
variations among the elements C, 0, Ne, Mg and Si were found to be 
relatively small — less than about ± 30 percent in the 5-15 
MeV/nucleon interval. Thus, propagation effects which would result in 
abundance time variations do not appear to account for the persistent 
difference — a factor of three to five — between the SEP and 
photospheric abundances of C and 0 relative to Mg and Si. 
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The velocity (or energy/nucleon) dependence of the SEP 
abundance ratios provides a more critical test for possible 
acccleration/propagation biases. This dependence may be derived from 
Figure 4.7 which shows, for each flare event, the measured 
energy/nucleon spectra for each of the elements C, 0, Ne, Mg. Si and 
Fe-fNi. Over the relatively narrow energy interval from 5 to 15 
MeV/nucIeon, the comparison of the spectral indices (slopes of the 
power law fits) suffices to accurately determine the energy/nucleon 
dependence of the relative abundances. In Figure 5.2, such 

comparisons are made for each flare event. Despite the variation of the 
spectral indices from flare to flare, the spectral indices of C, 0, Ne and 
Mg are roughly equal in a given flare event. However, the spectral index 
of Fe+Ni is often significantly different from those of C, 0, Ne and Mg. 
Tlie spectral indices of Si show some significant differences with 
those of C. 0, Ne and Mg, with a tendency in the direction of the Fe + Ni 
spectral indices. 

The relation between the SEP spectral inder and re'ative 
abundance measurements is explored in Figure 5.3. For every element 
pair in the set of elements (C, 0, Ne, Mg, Si and Fe+Ni) for which 
spectral indices were measured, the "relative abundance" is plotted 
versus the difference in spectral index measurements. Here, "relative 
abundance" refers to the relative abundance at 7.5 MeV/^nucleon, that 
is, the ratio of the "A" parameters of Figure 4.7. Also plotted are 
horizontal solid and dashed lines which indicate respectively the 
photosphcric relative abundance measurements and their estimated 
uncertainties. A number of interesting features are apparent in 
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Figure 5.2 

Spectral indices for the elements C, 0. Ne, Mg, Si and Fe-fNi in flare 
events i through vii. The (± la; error bars include the effect of counting 
statistics only. (7'and ore defined later in Section 5.3) 


SPECTRAL INDEX, y (ASSUMING dJ/dE oc E~ 
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Figure 5.3 

Relative SEP abundance plotted versus the difference in spectral index 
( 7 ,- 7 y) for the various element pairs among the elements C, 0, Ne, Mg. 
Si, and Fe+Ni. The (± la) error bars include the effect of particle 
counting statistics and, in the case of (7i”7y). a contribution due to 
possible systematic error (see Section 4.5.1), The horizontal solid and 
dashed lines indicate respectively the photospheric abundance results 
and their estimated uncertainties (from Meyer and Reeves 1977). 



oaAvoAo SEP EVENTS I through vii 
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Fi<?ure 5.3: 

t'l) The averaj{e value of 6y » 7x~7j is approximately xero for each 
element pair [x.y]. Thus, averaged over all seven flare events, the 
SEP relative abundances would be nearly independent of 
velocity (or energy/nucleon). This contrasts with the results of 
Crawford et. al. (1975) which suggested that in SEPs the heavier 
elements are always increasingly enhanced at lower energies. On the 
other hand, our result is consistent with the measurements of 
Mason, lirvestadt and Gloeckler (1979) which indicate that SEP 
composition at energies near 1 MeV/nucleon is nearly identical to 
that measured above 10 MeV/nucleon. The different conclusions of 
Crawford et. al. (1975) regarding the energy dependence of SEP 

composition may just be the result of a different flare sample. 
However, Mason, Hovestadt and Gloeckler (1979) have suggested that 
the low energy heavy element enhancements reported by Crawford et. 
al. (1975), which were obtained mainly in low altitude rocket flights, 
may ’ uv been caused by the earth’s magnetic field. (Also, see the 
revie. o*’ : vcent spacecraft observations by Mewaldt 1980). 

(2) There is essentially no correlation between 6y and relative 
abundance for any element pair. A correlation would be expected if 
the relative abundances were constant at some energy/nucleon Eq 
that was not too far removed from the observation range, 5-15 

MeV/nucleon, since in this case the relative abundance at 7.5 

MeV/nucleon would be proportional to (Eo/7.5)'’’'. For element pairs 

(like [Fe,X], where X = C, 0, Ne, Mg or Si) which show a large 
flare— lo-flare abundance variation, the lack of correlation between 
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<57 and relative abundance Indicates that there is no nearby 
energy/nucleon at which the relative abundances tend, to be constant 
from flare to flare. This also contrasts with the results of Crawford 
et. al. (1975) which suggested that the range of flare-to-flare variation 
of the SEP Fe/0 ratio was relatively small at energies above about 
15 MeV/nucleon, but larger at lower energies. 

(3) Among the different element pairs, the range of flare-to-flare 
variation of 6y and of the relative abundances are correlated, and tend 
to be larger for those element pairs in which the two elements have a 
larger difference of nuclear charge. For example, the spread of data 
points in both the vertical and horizontal directions increases as we 
move from the top left-most frame of Figure 5.3 downward, through 
the element pairs [C,0], [C,Ne] [C.Mg], [C,Si] and [ Fe]. 

(4) The range of flare-to-flare variations of 6j and of the relative 
abundances are not correlated with the size of the difference between 


the average 

SEP 

relative 

abundances 

and 

the photosphere 

abundances. 

For 

example. 

the element 

pair 

[Mg,0] has 6y 


approximately equal to zero in every flare event, and the Mg/0 ratio is 
nearly constant from flare to flare (only event ii shows a Mg/0 ratio 
that is significantly different from the Mg/0 ratio in the other events), 
yet the SEP Mg/0 values are a factor of four to five larger than the 
photospheric value. On the other hand, the element pair [Fe,Si] shows 
large flare-to-flare variations of both <57 and abundance ratio, while the 
average SEP Fe/Si ratio is nearly equal to the photospheric value. 
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In conclusion, the observations of this section suggest that the 
persistent differences which exist between the SEP and 
photospheric composition results are not due to systematic differences 
in SEP spectral index among the different elements, nor to systematic 
propagation effects which would result in time dependent abundance 
ratios. This suggests the other alternatives: (a) the composition of the 
flare acceleration site is significantly different from that measured 
for the photosphere and/or (b) there is a persistent acceleration bias 
which operates at energies outside our range of observation. 

5.3 The Selection of Four Preferred SEP Events 

For the purpose of studying solar composition, we concentrate on 
those flare events for which the measured SEP abundance ratios are 
approximately independent of energy/nulceon, and therefore have 
unique values which might equal those in the pre-accelerated plasma at 
the SEP acceleration site. 

For each flare event, a measure of the statistical significance 
of the differences among spectral indices of the various elements was 
calculated using: 

- T H(7j“7)Vaf 

^ j=i 

The elements C. 0. Ne. Mg, Si and Fe+Ni are indexed by j ; 7 j and Cj are 
the spectral ir.dex and its statistical error respectively (from Table 
4.12); ^ was chosen to minimize "^he results are listed in Figure 5.2. 
The four events with smaller values of (i. h. iv, and v) were chosen 
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for further study. The rejection of the other three events is based 
primarily on differences between the Fe+M spectral index and those 
of the other elements — differences which correspond to relatively 
large changes of relative abundance as a function of 

energy/nucleon. For example, the differences in spectral index of 
Fe+Ni and 0 (ype+Ni ” Ifo) of about +2 and -2, in events iii and vi 
respectively, correspond to factors of about ten («[15/'5]®) change in 
the (Fe+Ni)/0 ratio over the observation interval 5-15 MeV/nucleon. 

Note, in Figure 5.2, that the inclusion of flare event i in the set 
of preferred events is questionable because of the relatively low value of 
the Fe+Ni spectral index. Thus, in the next section we w'll consider the 
effect of event i on our conclusions. 

In a theoretical context, the flare events in which the SEP 
abundance ratios were approximately independent of energy/nucleon 
may be those in which the various nuclei, including Fe, were fully 
stripped of electrons' during acceleration. However, as discussed in 
Section 4.3, systematic time variations of the SEP (Fe+Ni)/0 ratio at 
5-15 MeV/nucleon, as well as direct charge state measurements at 
lower energies, indicate that the Fe nuclei are typically not ful'y 
stripped during their propagation through interplanetary space. For 
example, in flare event iv the measured energy/nucleon spectra are 
accurately the same among the different elements (including Fe+Ni and 
0). yet a systematic decrease of the (Fe+Ni)/0 ratio is observed early in 
the event. If we assume that, in general, the Fe nuclei are also not 
fully stripped during acceleration, then the above selection of flare 
events with SEP abundances independent of energy/nucleon may 
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correspond to the selection of events in which the important processes 
were primarily velocity dependent, rather than rigidity dependent. 

S.4 SystematicB of Flare-io-Flare SEP Abundance Variations for Four 
Preferred Events 

The systematics of the flare-to-flare composition variations are 
shown in Figure 5.4 by comparing the SEP composition results for each 
of the four selected flare events to the "four-flare average" 
composition. (The term "average" henceforth refers to the geometric or 
'log" average.) In each of the four flare events the deviations of the SEP 
abundances from their four-flare average values are approximately 
monotonic functions of nuclear charge Z in the range 6 ^ Z ^ 28. In 
particular, the flare-to-flare abundance variations of C, N and 0 are 
correlated, as are the abundance variations of Ca, Fe and Ni; while 
the abundance variations of C, N and 0 relative to Si are anticorrelated 
with those of Ca. Fe and Ni. An important exception to the 
approximate Z ordering of the abundance variations occurs for He whose 
abundance relative to Si is approxim_tely the same in all four events. 

The correlations of the SEP abundances seen in Figure 5,4 
suggest that the SEP compo.sition may be described by an average 
composition and a systematic deviation which varies in strength, but 
not character, from flare to flare. In particular, the ’Fe-rich" event ii 
does not stand out as a separate type of event, but rather may be the 
tail of a continuum of Z-dependent variability. 


The systematics of the flare-to-flare SEP composition variations 
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Figure S.4 

SEP flare-to-flare composition variations. For each of the flare events i, 
ii, iv, and v, the SEP abundances (from Tables 4.2 and 4.3 and normalized 
at silicon, ■ ) are divided by the four-flare average abundances (from 
Table 5.1) and are plotted versus nuclear charge Z. 



4 8 12 16 20 24 28 

ATOMIC NUMBER (Z) 
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seen in Figure 5.4 give practical significance to the four-fiare 
average composition, since we are able (in the next section) to compare 
the four-flare average abundances to those from other solar sources, 
knowing that the same comparisons, when made separately for the four 
flare events, would diflfer only by roughly monotonic functions of Z (for 
6S Z§ 20), Further, the average SEP composition is not sensitive to 
our choice of this particular set of four flare events. For example, 
tightening somewhat the meaning of "approximately independent of 
velocity" excludes event i, which has only a negligible effect on the 
average SEP composition. On the other hand, the average composition 
for the entire set of seven flare events is also not significantly 
different from the four-flare average. 

It is interesting that the selection of flare events with SEP 
abundance ratios that are approximately independent of 
energy/nucleon does not reduce the range of flare-to-flare variation of 
the Fe/Si ratio — the Fe/Si ratio ranges over a factor of nearly ten in 
both the full seven-flare and reduced four-flare sets. This is consistent 
with recent suggestions (Zwickl et. al. 1978; Briggs, Armstrong and 
Krimigis 1979) that the flare-to-flare variation of SEP elemental 
abundances may be largely due to variability of the composition of the 
pre-accelerated plasma at the SEP acceleration sites. 

5.5 Comparison of the Four-Flare Average SEP Composition to Other 
Solar Abundance Measurements 

The subject of solar composition was reviewed by Ross and Aller 
(1976) and updated recently by Aller (1980). The primary source of 
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abundance information is photospheric spectroscopy. Visible 
absorption line measurements, together with an atmospheric model and 
atomic transition probabilities, provide abundances for most of the 
elements — the major exceptions being Ne, Ar and the important 
constituent He. For most elements, photospheric spectroscopy is 
considered the most reliable source of solar composition 

information: (a) typical uncertainties in abundance measurements are 
estimated at less than about 30 percent, (b) the photosphere is 
thought to be well mixed with the outer convection zone of the sun by 
observed turbulent motions. 

Abundance information is also available for the solar corona. 
Extreme ultraviolet (EUV) and X-ray emission lines are used to 
obtain abundances of some elements, but with typically large 
uncertainties of a factor of two or three. A key question concerns the 
compositional homogeneity of the corona, that is. the role of 
turbulent mixing versus fractionation processes (see, e.g., Nakada 
1970, Wilhbroe 1976, and Mariska 1980). 

A third source of abundance information is solar wind 
measurements. Electrostatic deflection techniques provide abundances 
for H, He. Si and Fe (Same et. a.. 1975). The foil collection method 

gives He, Ne and Ar abundances (Geiss et. al. 1972). As with the corona, 
there is the possibility that the average solar wind composition may 
differ from that of the photosphere. The solar wind ratio ''aries 

with time by a factor of over ten. while the Si/H and Fe/H vary by 
somewhat smaller amounts (Same et. al. 1975). 
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Meteorite studies provide a fourth, Indirect, source of 
abundance information. For the non-volatile elements meteoritic 
abundances are in very good agreement with solar photospheric values, 
consistent with the idea that, for these elements, both the meteorites 
and solar photosphere reflect the composition of the primordial solar 
nebula (see Meyer 1978a, end Lambert and Luck 1970). Unfortunately, 
the most abundant solar constituents H, He, 0, and C are volatile and, 
as a result, are depleted in meteorites. 

In Figure 5.5 (also see Table 5.1) we compare the four-flare 
average SEP composition with results for the photosphere, corona and 
solar wind. We have taken the photospheric abundance data from 
Meyer and Reeves (1977) which is the most recent compilation which 
covers all the elements of interest here and incorporates error 
estimates. The photospheric abundance data of Meyer and Reeves 
(1977) arc in close agreement with the earlier compilation of Ross and 
Aller (1976) and with the recent photospheric abundance determinations 
by Lambert (1978) and Lambert and Luck (1978) (although Lambert 
and Luck estimate a significantly smaller uncertainty for the A1 
abundance). Figure 5.5a shows that relative to the photospheric 
abundance data the average SEP abundances of C, N, 0, and to a lesser 
extent S, are depleted with respect to the other elements in the range 
ll^ZgSB. The relative depletion of C, N and 0 in SEPs ic now seen 
to be a persistent effect, as reported in an earlier account of this work 
(Cook et. al. 1979) and by McGuire, von Rosenvinge and McDonald 
(1979). The depletion is also present in the earlier data of Teegarden, 
von Rosenvinge and MacDonald (1973), Crawford et. al. ( 1975 ), 
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l^igure 5.5 

Comparison of the four-flare average SEP abundances (normalized at 
silicon, ■ ) to abundances from (a) the photosphere (• , Meyer and 
Reeves 1977), (b) the corona (• , Meyer and Reeves 1977, C Parkinson 
1977, n Withbroe 1975, and A inferred from the N/0 measurement of 
McKenzie et. al. 1978 and Parkinson’s 0/Si ratio), and (c) the solar wind 
( # Boschler and Geiss 1976). The (± la) error bars include the 
estimated uncertainty in the photosphere, corona or solar wind 
abundances, and the uncertainty due to particle counting statistics in 
the average SEP abundances. The (± la) estimated uncertainty in the 
average SEP abundances due to flare -to -flare abundance variations is 
indicated for carbon and iron by brackets. 


SW/SEP CORONA/SEP PHOTO/ SEP 
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Notes to Table 5.1 

(» ) The geometric mean of the SEP abundances from flare events i, I?, iv 
and V. The (± lo) uncertainties include the effect of particle counting 
statistics, but not that of flare-to-flare abundance variations. 

( + ) Meyer and Reeves (1977). 

(§) (a) Meyer and Reeves (1977), (b) Parkinson (1977), (c) Withbroe 
(1975). 

( ^ ) Solar wind abundances relative to silicon were Inferred from: He/H = 
0.04±0,01. 0/H = (5±2 )x 10-4, Fe/H = (5±2)xl0“®. Si/H = 

(7.6± 3) X 10"^. He/Ne = 530±70, and Ne/Ar = 41± 10 (Boschler and 
Geiss 1976) 


( ^ ) Numbers in parentheses indicate factors of error. 
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Pellerin (1975), Webber et. al. (1975) and Nevatia, Durgaprasad and 
Biswas (1977). The only SEP composition measurements which do not 
show the depletion are the earliest nuclear emulsion results (e.g., 
Bertsch, Fichtel and Reames 1972) which, according to Webber et. al. 
(1975), could be in error as the result of poor charge resolution 
and/or detection efficiency. 

The comparison of coronal abundances (Figure 5.5b) from the 
compilation of Meyer and Reeves (1977) to the average SEP composition 
shows a relative depletion of SEP C, N and 0 which is less pronounced 
than that seen in Figure 5.5a. However, several recent coronal 
measurements show good agreement with the average SEP 
composition, most notably for 0 where large persistent depletions 
occur in SEPs relative to the photospheric data. Abundances of 0 
which are low relative to the photospheric results, but consistent 
with the SEP values, have been found in both EUV (Withbroe 1975; 
Flower and Nussbaumer 1975; Mariska 1980) and X-ray studies (Acton, 
Catura and Joki 1975; Parkinson 1977). 

Figure 5.5c shows that the solar wind and average SEP 
elem.3ntal abundances are in good agreement. Further, Figure 5.1 shows 
that the ranges of variation of the 0/Si and Fe/Si ratios are similar in 
the SEPs and solar wind, although some of the variation of the solar 
wind measurements may be due to experimental errors — see Bame et. 
al. (1975). (While the average elemental compositions measured for 
SEPs and solar wind are similar, there is a puzzling difference of nearly 
a factor of two between the SEP and solar wind measurements of the 
isotope ratio ^°Ne/^^Ne — Dietrich and Simpson 1979 and Mew'aldt et. 
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al. 1979.) 

In Figure 5.6 we focus on the 0/Si ratio where measurements are 
available from all four sources (SEPs, photosphere, solar wind and 
corona). The individual coronal measurements (EUV data, Mariska 
1980) are for a wide range of different coronal environments — quiet 
sun, coronal hole, active region and prominence. The 0/Si ratios 
measured in SEPs, solar wind and the corona have a comparable spread 
and, on the average, are all low relative to the photospheric value by 
slightly more than a factor of three. Although the 0/Si ratios 
observed in SEPs, solar wind and the corona may all be low for 
different reasons — e.g., possible preferential acceleration of SEPs 
and solar wind; and possible systeraatic errors in the solar wind (Same 
et. al. 1975) and coronal measurements (Meyer 1970b, Meyer and 
Nussbaumer 1979, Mariska 1980) — Figure 5.6 suggests a common 
cause. 

The overall similarity of the average SEP elemental 
composition and the elemental composition measured for the solar wind 
and corona, and in particular the evidence for a common persistent 
depletion of oxygen relative to the photospheric composition results, 
suggests that (a) the SEPs originate in the corona and (b) both the 
SEPs and solar wind sample a coronal composition which is persistently 
different from that measured for the photosphere. 

It is interesting that the Z-dependence of the ratios of 
pho tospheric abundances to four-flare average SEP abundances (seen 
in Figure 5.5a) is considerably different from the Z-dependence of the 


Figure 5.6 


Comparison of the 0/Si ratio measured in the photosphere, SEPs (this 
work), solar wind, and corona. For the SEP data, the filled circles refer to 
the four preferred SEP events (see Section 5,3). 
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SEP flarc-to-flare abundance variations (seen in Figure 5.4). The 
flare-to-flare SEP abundance variations were fairly smooth 
monotonic functions of nuclear charge Z, for 6^ 28. In contrast, 

the Z-dependence of the ratios of photospheric abundances to 
four-flare average SEP abundances shows a sharp break between 0 and 
Mg. and a feature at S — the Z-dependence is neither smooth nor 
inonotonic. 

While the ratios of average SEP to photosphere abundances are not 
ordered by Z, Figure 5.7 shows that they are roughly ordered by a 
different atomic parameter — the first ionization potential — as noted 
originally by Webber (1975). Restriccing attention to the portion of 
Figure 5.7 left of the dotted line, we see that the elements with first 
ionization potential less than 8 eV form a group in which there is 
rough agreement between the four-flare average SEP and 
photospheric abundances. However, the elements (C, N and 0) wdth 
first ionization potential above 1 1 eV are depleted with respect to the 
elements with ionization potential below B eV by a factor of about five. 
Sulfur, with an ionization potential near 10 eV, appears in the 
region of transition between the two groups of elements. 

The abundances of the other elements (Ne. Ar and He) with first 
ionization potential above 1 1 eV are not measured in the 
photosphere. However, on the right-hand side of Figure 5.7 we compare 
the average SEF’ abundances of these elements to the best available 
estimates (based primarily on observations of interstellar gas and 
hot stars) of their solar abundances from the compilation of 'Local 
Galactic" abundances by Meyer (1978a). The average SEP abundances 
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Figtire 5.7 

Ratios of the four-flare average SEP abundances to photospheric 
abundances (left of the dotted line) or to 'Local Galactic" abundances 
(right of the dotted line) plotted versus first ionization potential. All 
abundances are normalized to silicon (®). The "SEP uncertainty" 
is the (± Icr) uiu crtainty in the four-flare average abundance, including 
both the effects of particle counting statistics and flare-to-flare 


abundance variations. 


ERTAINTY 
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of Ne. Ar and He appear depleted similar to those of C. N and 0. 
However, since for Ne, Ar and He, we are not comparing the average SEP 
abundances to solar measurements, the right hand side of Figure 5.7 
should be viewed with caution. 

Figure 5.7 suggests that (a) first ionization potential or some 
other related atomic property (such as the cross section for ionization 
by electron impact) plays an important role in the SEP acceleration 
process, and/or (b) this same property is involved in the chemical 
differentiation of the corona from the photosphere. The first 
ionization potential is of major importance in the photosphere where 
normally temperatures are of the order of 10^ K. Here, elements 
v.ith ionization potential less than 7 eV are predominantly ionized 
while elements with ionization potential greater than 13 eV are mostly 
neutral (e.g., see Class 1951 and Gingerich, Noyes and Kalkofen 1971). 
Thus, neither possibility (a) nor (b) would be surprising. For example, 

(a) might occur if the SEP acceleration region is initially at the 
relatively low photospheric temperature. Then the electromagnetic 
acceleration of elements with high first ionization potential would be 
delayed until they are ionized and this delay might cause the depletion 
of elements with high first ionization potential seen in Figure 5.7. Case 

(b) — which is suggested by the similarity of SEP, solar wind and 

coronal composition measurements discussed above — could conceivably 
result irom the downward drift of neutral atoms in the 

photosphere— corona transition region. Neutral atoms experience the 
downward gravitational acceleration, but are not affected by the 
magnetic and electric forces which may support and accelerate the 
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charged component of the atmosphere. (The chemical difTerentiation of 
the corona has received some theoretical attention — see, e.g., Geiss, 
Hirt and Leutwyler 1970, Nakada 1970, and Shine, Gerola and Linsky 
1975. However, detailed predictions of coronal composition which could 
be compared to our SEP abundances are apparently not available,) 

Naturally, the correlation seen in Figure 5.7 must be viewed with 
caution since first ionization potential is related to the overall 
atomic structure. A search for other atomic parameters which 
may order the photosphere~SEP composition differences is under 
way —the results will be reported in a later paper. In particular, the 
search involves the parameters Z*/A and (Z*)^/A, where Z* — the 
average ionic charge of the nuclear species with atomic weight A — is 
taken from ionization equilibrium calculations (e.g. Jordan 1969) 
appropriate to coronal conditions and is a function of coronal 
temperature. The parameter Z*/A, the charge to mass ratio, may be 
important if the SEP-photosphere difference results from 
biased electromagnetic acceleration, while the parameter (Z*)^/'A 
enters if frictional drag due to Coulomb interactions is important. 

It is necessary to mention that the correlation seen in Figure 5.7 
depends on the correctness of the photospheric abundance 
determinations. Meteoritic abundances verify the photospheric results 
only for the non-volatile elements, which are exactly those in the group 
with first ionization potentials below 8 eV. 
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5.6 Summary 

The Low Energy Telescopes on the Voyager spacecraft have been 
used to measure the elemental composition (2 S Z S 28) and energy 
spectra (5 to 15 MeV/nucleon) of solar energetic particles (SEP) in 
seven large flare events. Aiming to understand the relationship 
between SEP elemental compc.<ition and the comp osition of the sun 
we have (a) discussed thf” Lime and velocity dependence of SEP 
composition in the individual fiore events, (b) select -J for further study 
four events with SEP abundances approximately independent of 
energy/nucleon, (c) examined the SEP flare-to-flare composition 
variations among the four selected flares and (d) compared the SEP 
composition for the selected events to solar abundance results for the 
photosphere, corona and solar wind. 

For the four selected events, the SEP composition results could 
be described by an average composition plus a systematic deviation 
about the average. In particular, for each of the four events, the ratios 
of the SEP abundances to the four-flare average SEP abundances were 
seen to be approximately monotonic .functions of nuclear charge Z in 
the range 6g 28. An important exception to this pattern of Z- 
dependent flare-to- flare abundance variation occurred for He, whose 
abundance relative to Si was nearly the same in all four events. 

The four-flare average SEP composition was found to be 
significantly different from the solar composition determined by 
photospheric spectroscopy. The ratios of the four-flare average SEP 
abundances to photospheric abundances are roughly ordered by first 


154 


ionization potential — the atomic parameter that determines which 
elemental species are ionized and which are neutral in the 
photosphere. Compared to photospheric abundance results, the 
elements with first ionization potentials above 11 eV (C. N and 0) are 
depleted in SEPs by a factor of about five relative to the elements with 
first ionization below 0 eV (Na, Mg, Al, Si, Ca, Cr, Fe and Ni). The 
abundances of the elements Ne, Ar, and He (which also have first 
ionization potential greater than 11 eV), are not measured in the 
photosphere. However, relative to 'Local Galactic" abundances of 
these elements, the four-flare average SEP abundances are depleted 
similar to those of C, N, and 0. 

For some elemental abundance ratios, the difference between SEP 
and photospheric results is persistent from flare-to-flare (even 
among all seven events) and is apparently not due to systematic 
differences in SEP spectral index between the elements, nor to 
propagation effects which would result in a time-dependent abundance 
ratio in individual flare events. A striking example occurs for the Mg/0 
ratio. The elements Mg and 0 have nearly equal spectral indices in each 
of the seven flare events; the Mg/0 ratio is approximately constant with 
time in each event and has nearly the same value from event to event. 
However, the SEP Mg/0 ratio is about a factor of four larger than the 
photospheric result, which has an estimated uncertainty of only 
about 25 percent. This suggests that (a) the composition of the flare 
acceleration sites are persistently different from that of the 
photosphere and/or (b) there is a persistent SEP acceleration bias which 
operates at energies outside our range of observation. 


The four-flare average SEP composition was found to be in good 
agreement with solar wind abundance results and with a number of 
recent coronal abundance measurements. This favors possibility (a) 
above. In particular, the evidence for a common depletion of oxygen in 
SEPs, the corona and the solar wind relative to the photosphere 
suggests (assuming the correctness of the photospheric results) that 
the SEP’s originate in the corona and that both the SEPs and solar 
wind sample a coronal composition which is significantly and 
persistently different from that of the photosphere. 
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Appendix A 
Data Proceaaing 

Data from the Voyager spacecraft are received at earth by the Deep 
Space Network and are processed at the Caltech Jet Propulsion 
Laboratory, which delivers "experiment" data (magnetic) tapes to the 
various Voyager science teams. The Cosmic Ray Subsystem (CRS) data 
are sent to the Goddard Space Flight Center and are reformatted to 
produce CRS "encyclopedia" data tapes, copies of whic.h are sent to 
Caltech. The tapes are logged and merged to produce a set of 
encyclopedia tapes containing time-ordered data. 

The data processing for this thesis is shown schematically in Figure 
A. 1. The encyclopedia tapes were processed using the VSTRIP program 
(Aufrance and Garrard 1978) to create the "strip" tapes VlSZOl (Voyager 
l) and V2SZ01 (Voyager 2). These tapes contain LET Z^ 3 and helium 
PHA events, together with complete CRS rate and status information, 
for each hour of data from September 1977 to May 1978. The "strip" 
tapes were processed by the program ELMER which produced (1) a 
printed summary of the data for each six hour period. Including 
preliminary PHA event counts for the most abundant elements in 
selected energy/nucleon bins and rate information needed to normalize 
the event counts to obtain preliminary flux measurements, (2) the 
same summary information on tape, and (3) a "event" tape containing 
only PHA events. 

An assortment of programs made use of the Zk 3 PHA events to 
complete the LET electronic calibration and deduce an adequate 
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generalized range-energy relation for the elements carbon through 
iron. These results and the Ll and L2 detector pathlength 
measurements were used by the program ZCAL to produce a "charge" 
tape which contains for each Z^3 PHA event; (1) a number indicating 
th 3 time of occurrence, (2) a number indicating which LET telescope 
the evert was from and whether or not the L3 detector was triggered, 
(3) the energy (in MeV) measured in each detector, and (4) the 
charge measurements Zl and, if L3 was triggered, Z2. The program 
HEMAbo produced a similar 'helium mass" tape which for every helium 
PHA event contains the above items (1) through (3) and the mass 
measurements Ml and. if L3 was triggered, M2. 

A variety of specialized computer programs were Uoed to 
investigate the data on the "charge" and "He mass" tapes, and to 
compute the flux and abundance measurements presented in Chapter 
4. In particular, for each flare event, the program HISTPL was used to 
produce Zl and (Zl+Z2)/2 charge histograms (e.g. Figures 3.9 and 
3.10) from which the 'PHA event" averaged abundances of Z^ 3 nuclei 
(Tables 4.2 and 4.3) were obtained. The program ELSORT selected 
and accumulated PHA events in three hour intervals for the fluence 
measurements of Tables 4.4 through 4.11. 
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Appendix B 
Problems 

B.l LET B of Voyager 1 

The response of this telescope was unusual in that (1) the He mass 
measurement, Ml, showed a relatively large energy and time 
dependence, (2) a relatively large "background" of events was present in 
the charge range near fluorine, and (3) the LI versus L2+L3 element 
tracks had a uncharacteristically blurred appearance. The problem may 
be electronic instability, but was not .investigated in detail: instead the 
dat *‘rom this telescope were excluded from use in the measurements 
presented in Chapter 4, Tables 4.1 through 4.12 

B.2 LET C oi Voyager 2 

On April 1, 1978 the LI detector of this telescope experienced a 
simultaneous shift in energy calibration and increase in crant rate, 
perhaps as the result of a light leak in the thin A1 entrance window. 
Subsequent data from this telescope were not used in the 
measurements of Chapter 4, Tables 4.1 through 4.12. 

B.3 Pulse Height Multiplication in the L!^ 35 /.un Detectors 

A background effect specific to the iron grou]: nuclei was 
identified in Sections 3.3 and 3.5 and necessitated the fur' ho- study and 
special handling which is discussed here. The affected PHA events were 
clearly seen in Figure 3 . 7 , a plot of the charge measurements Zl versus 
Z2, as a clump of points near (Z 1 =26,Z2-2B). Inspection of the LI, L2 
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and L3 pulse heights of the events in this clump revealed a very 
specific signature: 

(1 ) The events lie directly on the LI versus L3 track for iron, but 
lie just above the L2 versus L3 iron track, indicating that the Ll and L3 
pulse heights are normal, but that the L2 pulse height is abnormally 
large. 

(2) The events have relatively small L3 pulse heights (El 3 less than 
about 170 MeV), indicating that the nuclei responsible for the events just 
barely penetrated detector L2. 

The effect also occurs for a fraction of the iron nuclei which 
barely penetrate the Ll detectors. In this case, the affected events may 
be seen as the unusual clump of points near Z1=2B in Figure 3,6. 

The above signature suggests a pulse height multiplication effect in 
the LET 35 /xm surface barrier detectors that is very similar to one 
observed earlier in the response of such detectors to fission fragi ients. 
The effect (discussed by Walter 1969) occurs for nuclei which deposit a 
large ionization charge density in the depleted silicon near the gold 
electrode. The large ionization charge density is thought to induce 
tunneling of additional carriers (electrons) from the electrode 
through the thin oxide layer which separates the electrode from the 
depleted bulk silicon. The orientation of the LET 35 /xm detectors is 
consistent with this multiplication hypothesis: The gold electrodes face 
toward the L3 detector, such that the ionization charge density near 
a gold electrode is largest for nuclei which just barely penetrate the 
detector, as in (2) above. 
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To further study the problem, three detector PHA events with 
charge measurements Z1 > 24 were re-analyzed. For each event the Ll 
and L3 energy measurements (Elj and Els) were used to calculate a new 
charge measurement Z' and obtain an estimate E '12 of the energy 
deposited in detector L2 by solving; 

Tli “ R(Eli +E'l2+El3 , Z' , M) — R(E'i;2+El3 , Z' . M) 

Tuj = R(E'i^+El3 . Z' , M) - R(El3 . Z' . M) 

taking 

M = 2.132 • Z' . 

R is the generalized range-energy function discussed in Section 3.4. 
and Tl 2 are the mean pathlengths for the Ll and L2 detectors 
respectively. The results are illustrated in Figure B.l , a plot of El 2 (the 
measured L2 energy) versus E'l 2 (the L2 energy inferred from the Ll 
and L3 energy measurements). Normal events lie along the diagonal, 
while events with abnormally large L2 energy measurements lie above 
the diagonal. (Events below the diagonal may be due to L2 edge effects 
and other processes discussed in Section 3,5.2) The abnormal events 
occur predominantly at large values of E'l 2 which can only be obtained 
by nuclei which barely penetrate L2, consistent with observation (2) 
above. Figure B.2 shows two histograms of Els/E'ljj; one histogram for 
events which barely penetrate L2 (Eij < IVO MeV) and one for the 
remaining events (E^a > 170 MeV). These histograms ar:,' suinmed over 
all the LETs used on both Voyagers and indicate that (l) about 15 
percent of the events with E^j less than 170 MeV have L2 energy 
measurements which are abnormally large, by an average of about 
12 percent and (2) the L2 energy measurements are essentially normal 
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Fi(pire B. 1 

A scatter plot of the L2 detector energy measurement (E^a) versus the 
L2 energy deposition (E'i; 2 ) inferred from the Ll and L3 energy 
measurements, for three -detector PHA events withZl>24. 
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Figure B.2 

Two histograms of the ratio of the L2 detector energy measurement 
(Els) to the L2 energy deposition (E'ls) inferred from the Ll and L3 
energy measurements. The histograms are summed over 

three-detector PHA events with Zl> 24, from all the LETs used on both 
Voyager spacecraft. The histogram on the left includes only PHA 
events corrs.sponding to nuclei which just barely penetrate detector L2 
(Eu< 170 MeV), while the histogram on the right includes the remaining 
events (E>^> 170 Me\'). 
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for events with greater than 170 MeV. 

The special data analysis necessitated by the multiplication effect is 
discussed below. 

B.3.1 Three Parameter Analysia for Iron Nuclei 

The results for iron nuclei presented in Table 4.8 (b.7-15.0 
MeV/nucleon) are based on PHA events selected as follows: 

(1) Normal iron events (i.e. those with normal L2 pulse height) 
were selected using the charge consistency requirement and charge 
boundaries discussed in Section 3.5, and were counted if the total 
energy measurement, Elj+El 2 +El 3 . was appropriate to the 8.7-15.0 
MeV/nucleon interval. 

(2) Iron PHA events having abnormally large L2 pulse heights were 

identified as those events with charge measurements near 
(Z 1 =26,Z2=28). (Specifically, the requirement was (ZI-fZ2)/2 > 25 

and (Z1-Z2) < -0.85.) For each of these PHA events the L2 energy 

measurement was divided by 1.125 (the mean shift of the abnormal L2 
energy measurements: see Figure B.2) to obtain a corrected value E^. 
Then the event was counted if E^i +E i^ + El 3 was within the total energy 
interval corresponding to 8.7-15.0 MeV/nucleon. 

PHA events with abnormal L2 pulse heights constitute about 15 
percent of the iron nuclei counts presented in Table 4.2. Any 
additional systematic error (beyond that discussed in Section 4.4.1) in 
the iron nuclei counts due to the abnormal L2 pulse heights should be 


167 


small. For example, if no corrections at all are applied to the 
abnormal L2 energy measurements (as was the case for the iron 
abundances reported earlier in Cook ct. al. 1979 and Cook, Stone and 
Vogt 1900) , the resulting iron abundance measurements differ from 
those of Table 4.2 by less than 5 percent. 

B.3.2 Two Parameter Analysis for Iron Nuclei 

In the two parameter analysis of iron (for the fluence 
measurements of Tables 4.4 through 4.11) the multiplication effect in 
the LI detectors was taken into account as follows: 

(1) Iron PHA events with normal Ll pulse heights were 
identified as those events with charge measurement Zl in the range 
24.0 to 27.2 and were binned according to energy/nucleon as described 
in Section 3.4. 

(2) PHA events with charge Zl in the range 27.2 to 32.0 and L2+L3 
energy measurements of less than 170 MeV are primarily due to iron 
nuclei with abnormally large Ll pulse heights, rather than nickel or 
other nuclei. Therefore, the Ll energy measurements of these events 
were divided by 1.125 before energy/nucleon binning. The remaining 
events with Zl in the range 27.2 to 32.0 (i.e., those with El 2 +Etj > 170 
MeV) are mainly nickel nuclei and were included with no Ll energy 
adjustment. All PHA events were sorted into energy/nucleon bins as 
if they were due to iron nuclei, exactly as in (l) above. 

In the two parameter iron analysis, no special treatment was 
applied to account for multiplication effects in the L2 detectors. They 
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do not significantly afTect the Zl charge measurement and at most shift | 

a small fraction (about 5 percent) of the PHA events which should 
fall in the 7. 8-0. 7 MeV/nucleon interval into the 8.7-10.6 MeV/nucleon 
interval. 

The inclusion of nickel and other nuclei (27 § Z S 32) increases 
the fluence measurements over measurements of pure iron by roughly 
the combined abundance of these elements relative to iron, or about 10 
percent. However, the inclusion of these nuclei probably has only a 
negligible effect on the spectral index (7, see Section 4.5), since the 
energy spectra of the various elements of the iron group are likely to be 
similar. (Spectral index variations among the elements appear to be 
roughly ordered by nuclear charge Z. such that, In a given flare event, 
neighboring elements have similar spectral indices; see Section 5.2. 

Further, the .nickel to iron ratio [8.7-15 MeV/nucleon, Table 4.2] is 
nearly constant from flare to flare.) 

The multiplication effect in the Ll detectors is potentially an 
additional source of systematic error (beyond those discussed in 
Section 4.5) in the iron fluence measurements of Tables 4.4 through 
4.11 and the iron spectral index measurements of Table 4.12. Upper 
limits to these possible additional errors were obtained by 
recomputing the fluences and spectral indices without the Ll energy 
corrections discussed in (2) above. The changes in the iron fluence 
measurements were typically less than 10 percent for energy bins in 
the range 5. 0-8. 7 MeV/nucleon. and were zero for the other higher 
energy bins. The iron spectral indices changed by less than 5 
percent. Since the actual systematic error induced in the iron 
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spectral indices is probably much less than 5 percent (and therefore is 
small compared to statistical error and the possible systematic error 
from other sources discussed in Section 4.5) this error is negligible. 
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